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ABSTRACT 

We present the results of an extensive Hubble Space Telescope {HST) imaging study of ^100 Swift long- 
duration gamma-ray bursts (LGRBs) spanning 0.03 ^ z ^ 9.4 using relative astrometry from ground- and 
space-based afterglow observations to locate the bursts within their host galaxies. Using these data, we mea¬ 
sure the distribution of LGRB offsets from their host centers, as well as their relation to the underlying host 
light distribution. We find that the host-normalized offsets of LGRBs are more centrally concentrated than 
expected for an exponential disk profile, (R/Rh) =0.67, and in particular they are more concentrated than the 
underlying surface brightness profiles of their host galaxies. The distribution of offsets is inconsistent with the 
distribution for Type II supernovae (SNe) but consistent with the distribution for Type Ib/c SNe. The fractional 
fiux distribution, with a median value of 0.75, indicates that LGRBs prefer some of the brightest locations in 
their host galaxies but are not as strongly correlated as previous studies indicated. More importantly, we find a 
clear correlation between the offset and fractional fiux, where bursts at offsets R/Rh ^0.5 exclusively occur at 
fractional fluxes > 0.6 while bursts aX R/Rh ^ 0.5 uniformly trace the light of their hosts. This indicates that 
the spatial correlation of LGRB locations with bright star forming regions seen in the full sample is dominated 
by the contribution from bursts at small offset and that LGRBs in the outer parts of galaxies show no preference 
for unusually bright star forming regions. Finally, we find no evidence for evolution from z Z 1 to z 3 in 
the offset or fractional fiux distributions. We conclude that LGRBs strongly prefer the bright, inner regions of 
their hosts indicating that the star formation taking place there is more favorable for LGRB progenitor produc¬ 
tion. This indicates that another environmental factor beyond metallicity, such as binary interactions or IMF 
differences, may be operating in the central regions of LGRB hosts. 

Subject headings: gamma-ray burst: general 


1. INTRODUCTION 

Long-duration gamma-ray bursts (LGRBs) are the most en¬ 
ergetic explosions known in the universe, with a volumetric 
rate of < 1% of the core-collapse supernova rate. A wide 
range of studies of their accompanying afterglow emission 
and host galaxies have been used to shed light on the prop¬ 
erties of the bursts and their progenitors. In particular, the 
association with broad-lined Type Ic supernovae (Type Ic-BL 
SNe) and exclusive locations in star forming galaxies firmly 
established that LGRBs result from the deaths of massive stars 


are sub ject to a wide range of possible binary interaction sce- 
narios (|Podsiadlowski et al.||2004[ [van den Heuvel & Yoon| 
[2007] ). 

Since direct observations of LGRB progenitors are unlikely 
due to their low volumetric rate, insight into the nature of 
the progenitors has to rely on studies of the environments 
in which LGRBs occur. On a large scale, LGRB hosts have 
been found to generally be blue, compact, low luminosity, low 
mass, low met allicity galaxies with high specific star forma¬ 
tion rates (e.g.,|Le Floc’h et al.|2003 [Christensen et al.|2004 


established that LGRBs result from the deaths of massive stars , tion rates (e.g., | L e Jlo c^ h et alJ 2003 [ [ Christensen et al. | 2004 [ 

(e.g.JChristensen et al.|2004l|Woosley & Bloom|2006l|Wain-| [ Wainwright et al.|2UU7| | Savaglio et especiall^;^n 

wrigllt et alFo()7K Detailed observation! of the host galax- compared to core-collapse SN hosts ( jSvensson et al. |2010^. 

u oloA o These results have led to the notion of a preference for low 


[wright et al.[|2007 ). Detailed observations of the host galax 

ies have als o indicated a preference for low metallicity en 
vironments ([Stanek et al.|[2006t [Levesque et al.[[2010b [Gra 


ham & Fruchter||2013|), although some bursts occur m en 

vironments with solar metallicity ( [Levesque et al.[[2Q10a|b 
[Levesque] [2012] ). Similarly, afterglow observations estab 


fished that LGRBs are powere d by relativistic j ets with an 
energy scale of ^ 10^^ erg (e.g.. Frail et al.|20Ql| . These ob¬ 
servations support the idea that tne progenitors of LGRBs are 
rapidly-rotating massive stars that undergo core-collapse to 
form a hyper-accreting black h ole (collapsars: [Woosley|1993 ; 
[MacFadyen & Woosley||1999[ ). However, it is unclear at the 
present whether these massive stars are single or whether they 
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compared to core-collapse SN hosts ( jSvensson et al.[2010| ). 

These results have led to the notion of a preference for low 
metallicity progenitors. 

The sub-galactic environments of LGRBs can provide ad¬ 
ditional clues about the nature of the progenitors. To date, 
two major approaches have been employed in the literature: 
(i) measuring the offset of LGRBs relative to the centers of 
their hosts, and (ii) measuring the fractional brightness at the 
LGRB positions relative to the o verall distribution of light 
within their hosts. In this context, [Bloom et"ar| ( [20021 ) com¬ 
bined ground-based afterglow observations and Hubble Space 
Telescope {HST) host observations to measure the offsets of 
20 LGRBs, and found that they were consistent with being 
drawn from an exponential disk profile. This led them to con¬ 
clude that LGRBs were associated with star formation, and to 
rule out the compact object merger progenitor scenario. How¬ 
ever, the small sample size and the large uncertainties in many 
of the measured offse ts prevented a more detailed analysis. 
Fruchter et al.[ ( [2006| ), on the other hand, argued that the ir- 
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regular morphologies of at least some LGRB hosts necessitate 
the use of the fractional brightness technique (hereafter, frac¬ 
tional flux). Using this analysis for 30 LGRB hosts with HST 
observations, in comparison to 16 core-collapse SNe from the 
GOODS survey, these authors found that LGRBs are more 
highly concentrated in the brightest regions of their hosts. 
Svensson et al. ( |2Q10| ) expanded the sample and reached a 


similar conclusion, including that LGRBs occur in regions 
with a higher surface luminosity compared to core-collapse 
SNe. These results are in tension with the conclusion of 
Bloom et al.| p002| ) that LGRBs track an exponential light 


distribution. 

In addition to the tension between the results of previous 
studies, and the relatively small samples used in these stud¬ 
ies, it is important to note other potential challenges. First, 
the aforementioned studies were based on LGRB samples col¬ 
lected from multiple satellites using different trigger criteria, 
potentially introducing selection effects that are difficult to 
quantify. Second, the redshifts were mainly obtained from 
direct host galaxy spectroscopy (emission lines), leading to a 
bias toward lower redshifts and more luminous hosts; in con¬ 
junction with the small samples this also limited the ability to 
probe any redshift evolution in the LGRB locations. Finally, 
these studies did not include “dark” bursts (i.e., bursts that 
suffer from large rest-frame extinction), leading to a potential 
bias against dusty environments; at least on the global scale 
the hosts of dark bursts appear to be mor e luminous and mas 
sive than the gene ral LGRB host sample ( [Kruhler et al.|201 1 


[Perley et al.|2013| ). These shortcomings can now be overcome 
with the much larger and uniform sample of Swift LGRBs. 

Beyond the direct study of LGRB progenitors, the frac¬ 
tional flux distribution of LGRBs has also been used as a 


point o f comparison with core-collapse SNe. |Fruchter et al. 
( 2006| ) found that the core-collapse S Ne from GOODS were 
consistent with a uniform distribution. [Svensson et al. ( |2010j ) 
reached a similar conclusion using a larger sample of GOODS 
and PANS core-co llapse SNe. For more local SN samples, 
[Kelly et al.j (|2008|) measured the fractional fiux distributions 
for Type 11 and Ib/c SNe and found that Type Ic SNe are con¬ 
sistent with being drawn from the same fractional fiux dis¬ 
tribution as LGRBs, whereas Type II and Ib SNe uniformly 
track the light of their hosts. They suggest that LGRBs and 
Type Ic SNe share a common progenitor, and that a factor 
such as metallicity may determine whether core-collapse re¬ 
sults in a LGRB with an associated Ic-BL SNe or a normal 
Type Ic SNe with no associated LGRB. However, recent work 
has found that LGRBs and Ic-BL SNe occur in host galaxies 
with high star formation rate density and that th is preference 
cannot be due to metallicity ( [Kelly et ar][2Q14[ ). Finally, the 
LGRB fractional fiux distribution has also been compared to 
that for super-luminous supernovae yielding a result that is 
suggestive of agreement ( [Lunnan et ai.|2015] ). With the grow¬ 
ing samples of SNe, comparative studies of SN and LGRB 
environments have become limited by the small sample of 


LGRBs from the studies carried out a decade ago (Bloom et 
[ar][2Q02[ [Fruchter et al.[[2Q06l ). This situation can now be 
remedied with the much larger sample of Swift LGRBs. 

Here we present the first uniform analysis of HST follow-up 
observations of the host galaxies of about 100 Swift LGRBs 
collected over the last decade. The goal of this analysis is 
twofold. First, to investigate the offset and fractional fiux 
distributions for a much larger and more uniform sample of 
LGRBs than previously possible. Second, to provide a much 
larger comparison sample for future studies of the locations of 


other astrophysical transients, such as various SN types. The 
structure of our paper is as follows. In Section 2 we present 
the details of the host galaxy and afterglow observations, data 
analysis, and astrometric matching. In Section 3 we discuss 
our offset, fractional fiux, and host assignment methodolo¬ 
gies, and present our measurements. In Section 4 we present 
our resulting offset and fractional fiux distributions, and in 
Section 5 we discuss potential trends with redshift, the rela¬ 
tionship between fractional fiux and offset, and implications 
for LGRB progenitors. We conclude with a summary in Sec¬ 
tion 6. 

In this paper we use Hq = 67 km s~^ Mpc~\ = 0.32, 
and Oa = 0.68 ([Planck Collaboratio n et al.|2014[ ), as well as 
AB magni tudes IQke & Gunn[|1983|), c orrected for Galactic 
extinction ( [Schlafly & Finkbeiner|201l] ). 

2. OBSERVATIONS AND ASTROMETRY 
2.1. HST Observations and Reduction 

To locate the LGRBs within their host galaxies we utilize 
high resolution HST images. We primarily use Wide Field 
Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) 
data, but also Wide Field Planetary Camera 2 (WFPC2) and 
Near Infrared Camera and Multi Object Spectrometer (NIC- 
MOS) data when the former are not available. Our sample is 
composed of 100 LGRBs observed over the last decade (Au¬ 
gust 2004 to July 2013) under multiple programs (see Table H 
for program ID numbers). We perform the first uniform anaL 
ysis of these data and utilize them for the purpose of inves¬ 
tigating the locations of LGRBs within their hosts and their 
relation to the underlying host light distribution. This sample 
has the significant advantage of largely consisting of bursts 
discovered by the Swift satellite and thus does not suffer from 
any potential biases that may be associated with combining 
data from multiple satellites as was done in previous stud¬ 
ies. About 86% of the bursts in the sample have measured 
redshifts (0.03 < z < 9.4) and the HST imaging for these 
bursts span rest-frame UV and optical wavelengths (Figure 
[TJ. The wide redshift range allows us to investigate trends as 
a function of cosmic time. Unlike previous studies, the sam¬ 
ple also includes a substantial number of dark LGRBs, which 
either lack optical afterglows or are much fainter than ex¬ 
pected (GRBs 051022, 060719, 060923A, 061222A, 070306, 
070521, 070802, 080207, 080325, 080607, 081109, 081221, 
090404, 090407, 090417B, 090709A, 100413A, 100615A, 
and 110709B). Some are detected in the NIR while others are 
exclusively detected in the radio and X-ray bands. 

We retrieved the HST observations from the MAST archive; 
for the ACS images we obtained charge transfer efficiency 
(GTE) corrected images, while for WFC3/UVIS data we used 
software from STScI to apply GTE correc tions. We re¬ 
duced the data using the astrodrizzle ( [Gonzaga et al. 
20121 ) task in the STSDAS IRAF package utilizing recom¬ 
mended parameter settings for each instrument. By combin¬ 
ing dithered exposures, this task enables the reconstruction 
of a higher resolution image than sampled by the instrumen¬ 
tal point spread function. The task also applies distortion 
corrections to the images, which are critical for precise as¬ 
trometric alignment. We use f inal_pixf rac = 0.8 and 
f inal_scale = 0.065", 0.02", 0.03", 0.05", and 0.15" per 
pixel for WFC3/IR, WFC3/UVIS, ACS, WFPC2, and NIC- 
MOS, respectively. In Table we list for each LGRB the 
program ID, instrument, filter, and exposure time of the final 
drizzled image. In cases where there are multiple epochs of 
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Figure 1. Rest-frame wavelength versus redshift of the observations used to measure offsets and fractional flux for each LGRB. Only bursts with measured 
redshifts are included in this plot. Also shown are projected histograms of the redshift distribution and the distribution of rest-frame wavelengths probed by the 
HST observations. 


imaging we list only the final epoch, unless an earlier epoch 
was chosen as the best host image. In Figure we show the 
final drizzled image for each burst, with the location of the af¬ 
terglow shown as determined by relative or absolute astrome¬ 
try (Sections 2.3 and 2.4). 

2.2. Afterglow Imaging 

To precisely locate the LGRBs within their host galaxies 
requires the detection of an afterglow. About 85% of bursts 
in the sample have detected optical or NIR afterglows. We 
utilize the deepest, highest resolution afterglow images avail¬ 
able to best match the quality of HST imaging. We primarily 
rely on publicly available afterglow images from the 8-m Very 
Large Telescope (VLT) and 8-m Gemini North and South tele¬ 
scopes. In 13 cases, an afterglow image is available from 
HST. We also use images from the Swift UV/Optical Tele¬ 
scope (UVOT), 6.5-m Magellan telescopes, Palomar 60-inch 
Telescope, Swope 40-inch Telescope, and the 8.2-m Subaru 
Telescope. 

We use standard IRAF packages to analyze the afterglow 
data. We use the task uvotimsum in the HEASoft pack¬ 
age to co-add SwiftlMWOT images, combining multiple fil¬ 
ters if necessary to obtain a high signal-to-noise ratio (S/N) 
afterglow detection. For HST afterglow images we use 
astrodrizzle as described in Section 2.1. For GRB 


070306 we use the eclipse pipeline to reduce the NIR data 
from VLT/ISAAC. For GRBs 080325 and 090709A we use 
the MCSRED IRAF package to reduce the NIR data from 
Subaru/MOIRCS. For the remaining bursts without optical 
or NIR afterglows, we rely on radio or X-ray detections. 
We obtain the reported radio positions from the GCN cir¬ 
cular archive. When an X-ray afterglow image is available 
from Chandra, we retrieve the image from the Chandra Data 
Archive. For each burst, we list in Table the telescope and 
filter utilized to obtain the afterglow image that was used for 
astrometry. 

2.3. Relative Astrometry 

To precisely locate the afterglows on the HST images we 
use relative astrometry. We first identify common point¬ 
like sources betwe en the afterglow and H ST images, us¬ 
ing SExtractor ( Bertin & Arnouts||1996| ) to measure the 
source positions. We then match the coordinates with the 
IRAF task ccmap, using a second or third order polynomial 
to compute the astrometric solution. Using the afterglow im¬ 
ages from the telescopes listed in Section 2.2 we generally 
obtain the following astrometric tie uncertainties, dtie (Icr): 

• 6 mas using an HST detection of an optical or NIR af¬ 
terglow 
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(a) Panel One 


Figure 2. Hubble Space Telescope drizzled images of 100 LGRB host galaxies studied in this work. The circles (3cr) indicate the location of each GRB. North is 
up and East is to the left. Red, blue, green, cyan, and yellow circles indicate that afterglow positions were determined from optical/NIR, SwiftIVNOT, Chandra, 
radio, and XRT images, respectively. Physical scales are given when the LGRB redshift is known. 
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Figure 2. Continued 
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Figure 2. Continued 
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Figure 2. Continued 
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Figure 2. Continued 
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• 30 mas using an optical or NIR afterglow detection 
from ground-based telescopes 

• 100 mas using a UV or optical afterglow detection from 
Swiftl\]NOT 

• 300 mas using an X-ray afterglow detection from Chan¬ 
dra 

In a few cases, the lack of sufficient common sources for 
a robust astrometric tie necessitates the use of an intermedi¬ 
ate image to locate the afterglow position on the HST image. 
This is the case for bursts detected with Chandra. We use the 
CIAO task wavdetect to measure the positions of X-ray 
sources. In these cases, the uncertainty associated with the 
match between the Chandra and intermediate ground-based 
images dominates the total astrometric tie uncertainty. Table 
|2] lists the number of common sources and astrometric tie un¬ 
certainty for each burst. Bursts with no listed astrometric tie 
uncertainty are those for which the best afterglow position is 
from XRT or for which we were not able to obtain an after¬ 
glow image. GRB 060418 is a unique case where the after¬ 
glow and host are well-detected in a single HST image and are 
sufficiently offset from each other such that their positions are 
not biased by light contribution from each other (later epoch 
imaging revealed no host emission at the afterglow position). 
Therefore astrometric matching is not necessary for this burst. 


of the light distribution) determined by sequentially increas¬ 
ing the threshold for detection inSExtractor. The statisti¬ 
cal uncertainty on the centroid determined by SExtractor 
does not take into account the larger systematic uncertainties 
caused by the fact that LGRB hosts tend to have irregular mor¬ 
phologies. Instead, the uncertainty on the host center (cThost) is 
estimated by taking the quadrature sum of the standard devi¬ 
ation in both directions of the positions measured by multiple 
runs of SExtractor using varying thresholds for detection 
from 3(7 to the highest threshold that yields a detection. 

Using the afterglow position and host center we calculate 
the projected angular offset {R) between the afterglow posi¬ 
tion and host center. The total uncertainty on the angular off¬ 
set {(Tr) is taken to be the quadrature sum of (Jtie, (Tag. and 
cohost- We also calculate the projected physical offsets using 
the measured redshift for each burst. For bursts with no mea¬ 
sured redshift we use the fact that at z > 0.5, typical of most 
LGRBs, the angular diameter distance is approximately con¬ 
stant so that we can use a conversion factor of about 8 kpc/". 
The projected angular and physical offsets and corresponding 
uncertainties are listed in Table [J] Bursts for which no offsets 
are listed indicate either bursts wi th w ell-localized afterglows 
but no detected host (see Section |3.3| for a discussion of dis¬ 
tinguishing between no host detection and a large offset) or 
bursts for which we only have an XRT position or were un¬ 
able to access the afterglow image. The uncertainties on the 
XRT positions are typically > 1" in radius, too large to enable 
a reliable host association. 


2.4. Absolute Astrometry 

In cases where the afterglow is detected only in the ra¬ 
dio band, or when the X-ray observations do not contain any 
common sources for relative astrometry, absolute astrometry 
is performed by matching the late-time HST images to the 
2MASS astrometric system. This is mainly the case for dark 
GRBs. For the two GRBs where a radio position was used 
(090404 and 100413A), the uncertainty in the afterglow po¬ 
sition with respect to the HST image is dominated by the as¬ 
trometric accuracy of those radio positions, ^ 0.5". Similarly, 
when using absolute astrometry to locate the X-ray afterglow 
of GRB 090407, the error is dominated by the astrometric ac¬ 
curacy of Chandra, ^ 0.4". 


3. MEASUREMENT METHODOLOGY 
3.1. Offset Measurements 

Once the afterglow and HST images are registered, we iden¬ 
tify the host galaxy in the HST image. First we measure the 
afterglow centroid using SExtractor to determine the af¬ 
terglow position relative to the HST image. The afterglow 


positional uncertainty ((Jag) is estimated as 


FWHM 


In some 


2(S/N) • 

cases, it is apparent in the afterglow image that there is an un¬ 
derlying extended object, presumably the host galaxy. To reli¬ 
ably measure the centroid of the afterglow we p erform imag e 
subtraction, using the ISIS software package ( Alard|[2Q00] ), 
relative to a late-time image of the host ideally from the same 
telescope and in the same filter. Once the afterglow position 
is determined, we u se a uniform host assignment procedure 
discussed in Section l33l Here we continue our discussion of 
the offset measurement procedure assuming a host galaxy has 
been identified. 

We measure the host galaxy center and corresponding un¬ 
certainty using SExtractor. We use a uniform definition 
of the host center as the brightest region of the galaxy (or peak 


3.2. Host Galaxy Sizes, Normalized Offsets, and Magnitudes 

As our interest lies in the positions of LGRBs relative to 
the light of their host galaxies, we normalize the offsets by 
the sizes of the hosts. Here we use the half-light radius (Rh) 
as measured by SExtractor (with ELUX_RADIUS = 0.5). 
We also measure the 80% light radius (Rgo) for comparison 
with previous studies of LGRBs and SNe. In Table we 
list the half-light radii and host-normalized offsets (Rnorm) for 
each burst. 

We also measure the apparent magnitudes of the host galax¬ 
ies, and corresponding uncertainties, using the MAG_AUTO 
estimate from SExtractor which uses Kron apertures. We 
obtain the magnitude zeropoints from STScI where tabulated, 
or from the photometry keywords in the HST ima ge headers. 
The magnitudes corrected for Galactic extinction ([Schlafly & 


Finkbeine'r|2011| ) and corresponding uncertainties are listed 
in Table tor each burst with an associated host galaxy. 

3.3. Host Galaxy Assignment 

To assign host galaxies to each LGRB we calculate the 
probability of chance coincidence, for several candidate ex¬ 
tended objects surrounding the afterglow position, given the 
observed number density of field galaxies from deep optical 


and NIR surveys (jHogg et al^|1997t |Beckwith et al.||2006 
Metcalfe et al.|2006| ). Nominally we assign the candidate with 
)robabiliL 


the lowest probability of chance coincidence as the host. For 
candidates at small offsets from the afterglow position, such 
that the afterglow position is essentially coincident with the 
candidate, this candidate will in most cases have the lowest 
probability of chance coincidence. Bursts for which there is 
no coincident candidate require further investigation to distin¬ 
guish between the scenario where the host was not detected 
and the scenario where the burst occurred at a large offset 
from a neighboring galaxy. We discuss these cases in detail 
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in Sec tion 3.3. 1| Following the methodology of Bloom et al. 
( 2002| ) we determine the probability of chance coincidence 
(Fee) using the following equation: 


1-e 


-T^K a{<m) 


( 1 ) 


where cr(< m) is the observed number density of galaxies 
brighter than magnitude m, and Re is the effective radius. 
For optical observation s we use the ga laxy number densities 
given by Equation 1 of Berger ( |2010| ), whereas for NIR ob- 
servations we use the //-band number counts measured by 
[Metcalfe et al.| ( |2006) ). As discussed by jBloom et al.| ( |2002| ), 
the appropriate value of Re depends on the offset, size of the 
gal axy, and unce r tainty on the afterglow position. Follow¬ 
ing [Bloom et al. ( 2002[) we take Re to be the maximum of 

(3 , y -h 4F^). In other words, if the localization 


is precise, such that the first term is much smaller than the 
offset and galaxy size, then the second term dominates. On 
the other hand, if the GRB localization is poor the positional 
uncertainty dominates over the offset and galaxy size, and its 
radius becomes the appropriate value of the effective radius 
for determining Fee- In the majority of cases 85%) the lo¬ 
calizations are sufficiently precise such that the s econ d term 
dominates over the first. As discussed in Section [2.3[ dtie is 
typically ^100 mas for bursts where we use afterglow detec¬ 
tions from Swiftf\]yOT or Chandra. In these cases, the GRB 
positional uncertainty can dominate over the offset and galaxy 
size. 

We list the values of Fee in Table For the bursts with co¬ 
incident extended objects. Fee tends to fall in the range 10“^ 


to 5 X 10“^. Therefore we do not expect significant contami¬ 
nation of our host galaxy sample by unrelated galaxies along 
the line of sight. 


3.3.1. Bursts With No Coincident Host Candidate 


The subset of bursts that lack a coincident extended object 
must be examined further to determine if the object with the 
lowest Fee can realistically be considered the host galaxy. Re¬ 
lying simply on Pcc may result in the inclusion of unrealisti¬ 
cally luminous objects in the sample. We identify as possible 
host candidates the nearest extended objects with Fee ^ O-l- 
Given the measured redshift of the burst, we can determine 
whether the luminosity of the host candidate is consistent with 
the range of LGRB host luminosities observed at that redshift. 
It may be more likely that the true host lies beyond the magni¬ 
tude limit of our observations and was not detected. However, 
it is important to avoid biasing our results by eliminating po¬ 
tential large offsets. Here we discuss our assessment of bursts 
whose host candidates satisfy Pcc ^01 but are offset by sig¬ 
nificantly more than the half-light radius of the host candidate. 

To carry out the assessment of non-coincident host candi¬ 
dates with Fee ^0-1 we compare their luminosities to the lu¬ 
minosities of the coincident host galaxies. We convert our 
measured apparent magnitudes to rest-frame luminosities in 
units of L*. Given the redshift, we determine the absolute 
magnitude and using measured values of L* from galaxy sur¬ 


veys ([Faber et al.[[2007[ [Marchesini et al.[[2007[ [Bouwens et 

al .]2015[ ) we calculate L/L* for each host galaxy. For bursts 


with z ^ 1.7 we use the measured val ue of L* in the F-b and 
for the blue galaxy sample analyzed by [Faber et al.[ ( [2007 ). In 
the range 1.7 < z < 3.2 we use L* values in the F-band from 
Marchesini et al.[ ( |2007[ ) using their U-V < 0.25 mag galaxy 
subsample. For bursts with z > 3.2 we use L* corresponding 


to UV bands as measur ed by surveys of Lyman-break galaxies 
( [Bouwens et al.|20f5] ). These were chosen so that the wave¬ 
length regime in which the L* values are appropriate corre¬ 
spond roughly to the rest-frame wavelengths probed by our 
observations in these redshift bins. We do not perform a K- 
correction to a common band. 

We plot the luminosities as a function of redshift for the full 
sample, including coincident and non-coincident candidates, 
in the left panel of Figure The right panel shows the lumi¬ 
nosities as a function of host-normalized offset. The luminosi¬ 
ties of the robust host galaxy associations span ^ (0.01 -2)L*, 
whereas some host candidates at relatively large offset and 
high Fee have luminosities ranging from ^ (3-100)L*, in¬ 
consistent with the distribution for bursts with coincident host 
galaxies and indicating that they are unrelated galaxies. The 
luminosity of the host candidate for GRB 060526 is ^ lOOL* 
assuming it is at the redshift of the GRB, much larger than 
expected for a GRB host. We therefore reject this associa¬ 
tion. We also reject the candidates for GRBs 060927, 070125, 
071031, and 091208B because their luminosities are unrea¬ 
sonably high given the observed distribution for securely as¬ 
signed hosts. In the case of GRB 080928, there are two can¬ 
didates with Fee ^0-1 where the object with Fee ~ 0.023 has 
a luminosity of ~ 4L* and the object with Fee ~ 0.035 has a 
luminosity of ^ L*. We therefore accept the lower luminosity 
galaxy as the host. We also accept as the host the candidate 
for GRB 081008 due to its reasonable luminosity of ^ 0.6L*. 
The left panel of Figurej^also shows upper limits on the lumi¬ 
nosities of the hosts where we are confident the host was not 
detected due to the lack of any nearby host candidates with Fee 
<0.1. Upper limits for the cases with rejected host candidates 
discussed above are also shown. 

Furthermore, we note that there is no trend of host lumi¬ 
nosity versus host-normalized offset. There is no reason to 
suspect from predictions that bursts at large or small offsets 
occur only in bright or faint galaxies. As shown in Figure 
the inclusion of the rejected host candidates creates an artim 
cial trend because the only galaxies at higher normalized off¬ 
set that have Fee ^0-1 are bright galaxies. Therefore we are 
confident these rejected candidates are unrelated galaxies and 
the true host was not detected. Indeed, the resulting upper lim¬ 
its on the luminosities of these hosts are not surprising given 
the upper limits we calculate for unambiguous non-detected 
hosts. Upon completing our host assignment procedure we 
obtain a sample size of 71 LGRBs for which we make offset 
and fractional fiux measurements. 

3.4. Fractional Flux Measurements 

Studying the locations of LGRBs within their hosts can be 
expanded beyond an analysis of the offset distribution. In 
many cases where the host’s morphology is irregular, the lo¬ 
cation of the host center is unclear and the offset from an 
empirically defined center may not be meaningful. Further¬ 
more, the host galaxy size defined by the half-light radius does 
not capture the complex morphology of many of the hosts. 
While measuring offsets is useful for comparing the locations 
of LGRBs with predicted profiles of star formation, we would 
like to go one step further and ask whether or not LGRBs are 
spatially coincident with sites of massive star formation. This 
is possible using a method, less sensitive to galaxy morphol¬ 
ogy than offsets, of comparing the brightness at t he LGRB 
location re lative to the entire host light distribution ( [Fruchter[ 
[et al.[2006| ). We calculate the fraction of the total flux from the 
host galaxy that is contained in pixels fainter than or equal to 
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Figure 3. Left: Host galaxy rest-frame luminosity in units of L* as a function of redshift for assigned host galaxies (red), rejected host candidates (green), and 
non-detected hosts (blue). The data point area is proportional to Pec where a smaller point mean s a sma ller Pec- The red dotted line indicates a luminosity of L*, 
and the dotted cyan line connects the two possible host candidates for GRB 080928 (see Section|3.3.1). We do not include bursts with unknown redshift. Right: 
Host galaxy luminosity versus host-normalized offset. As expected this is a scatter plot for the robust host associations. The correlated tail caused by the rejected 
candidates is an indication that they are unrelated galaxies. 


the flux at the LGRB location. The resulting "fractional flux" 
is a statistic that measures the brightness of the burst site com¬ 
pared to the entire galaxy. A fractional flux value of 1 means 
that the burst occurred on the brightest pixel in the galaxy. 

We measure the fractional flux for each LGRB using the 
following procedure. First we determine the flux value at the 
LGRB position. When the la error circle deflned by the un¬ 
certainty on the LGRB position as determined by the quadra¬ 
ture sum of atie and aAc spans less than 1 pixel, we simply 
take the flux value of that pixel. When the error circle encom¬ 
passes more than one pixel but is smaller than the PSF, we 
measure an average flux within the la error circle weighted 
by the fractional area encompassed by a given pixel. These 
scenarios are relevant for most bursts where we use afterglows 
from ground-based optical/NIR telescopes and HST (see Sec¬ 
tion 2). Using SExtractor we extract the host galaxy pix¬ 
els using a threshold of la above the sky background and then 
calculate the fraction of the total flux in pixels fainter than or 
equal to the average flux of the LGRB error region (the frac¬ 
tional flux). 

When the error circle is larger than the PSF, the galaxy 
brightness as a function of position could vary signiflcantly 
so that the uncertainty in the position of the LGRB may have 
a large impact on the fractional flux. This scenario is of¬ 
ten encountered when using optical afterglow detections from 
Swiftl\]NOT, X-ray afterglows from Chandra, or radio after¬ 
glow detections, where the error regions can encompass tens 
to hundreds of pixels. In some cases, these pixels include both 
the faintest and brightest pixels of the host galaxy. In this 
regime, it is unclear whether the fractional flux determined 


as described above remains a meaningful quantity. Here we 
show that for bursts with error circles larger than the PSF, 
a different weighting is required to properly assess the frac¬ 
tional flux in these cases. 

To assess this, we employ a different procedure for deter¬ 
mining the fractional flux, taking into account the Gaussian 
error distribution associated with each burst. As before we 
extract the galaxy pixels using SExtractor. Each pixel in 
the image has an associated probability for hosting the burst, 
governed by the 2D Gaussian distribution that takes into ac¬ 
count the astrometric tie and afterglow positional uncertain¬ 
ties. For all pixels encompassed by the error circle, a corre¬ 
sponding fractional flux value can be calculated. Pixels not 
associated with the galaxy will have fractional flux values of 
zero. The posterior probability distribution for the fractional 
flux can then be constructed by weighting each value by its 
corresponding Gaussian probability. The fractional flux and 
corresponding uncertainty are estimated as the mean and stan¬ 
dard deviation of this distribution. 

By comparing the fractional flux values for bursts with large 
error regions with the well-localized bursts, we And that there 
exists a considerable bias towards lower fractional flux val¬ 
ues in the former. We hypothesize that not only is the abso¬ 
lute size of the error circle important but also the error circle 
size relative to the galaxy size. If the galaxy is much larger 
than the error circle, even a large error circle in an absolute 
sense can yield a well-constrained fractional flux value. If a 
relatively well-localized burst occurs in a small galaxy, the 
fractional flux value may not be well-constrained. Figure 
illustrates this bias by plotting fractional flux versus the ratio 
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Figure 4. Left: Fractional flux versus the ratio of error circle area to galaxy area. Green points with error bars represent bursts with error circles larger than the 
HST PSF. Right: The fractional flux distribution binned into two bins separated at an area ratio of 0.1. In both pl ots it is apparent that there is a bias to lower 
fractional flux values when the error circle covers a substantial fraction of the galaxy area. In the analysis in Section [4!2] we only use bursts with a ratio of < 0.1 


of error circle area to galaxy area (left panel). The right panel 
of Figure]^ shows the cumulative distributions for two subsets 
of the data divided at an area ratio of 0.1. Of the bursts with 
area ratios ^0.1, there are none with fractional flux values 
greater than ^0.6, while bursts with an area ratio of < 0.1 
span the full range of fractional flux values. This bias can be 
understood by noting that as the error region increases relative 
to the galaxy size, there is greater probability that the burst 
occurred on background sky pixels. In general we And that 
for error regions large in an absolute or relative sense there 
is a bias towards lower fractional flux values. To avoid bias¬ 
ing our results we remove from consideration the bursts with 
area ratios of > 0.1, resulting in a sample size of 50 LGRBs. 
By removing these bursts we are not concerned about select¬ 
ing against bursts with intrinsically low fractional flux values. 
Here we are removing a bias associated with the measurement 
of fractional flux when burst localization is poor. 

In Table [3] we list the fractional flux values and area ratios 
for each LGRB. Bursts with no fractional flux value are ei¬ 
ther bursts for which a host was not detected or a sufficient 
afterglow position was not available. 

4. RESULTS 

We have compiled the largest sample of LGRB offsets and 
fractional flux measurements presented to date. Our aim is 
to use them to understand where within their host galaxies 
LGRBs occur and what the results from a large sample indi¬ 
cate about the progenitors. Here we present the results of our 
offset and fractional flux distribution measurements summa¬ 
rized in Tables [2] and O 

4.1 . Offset Distribution 

In this section we present the measured physical and host- 
normalized offset distributions. We also present the results of 


an analysis of the offset uncertainties to understand how they 
affect the offset distributions. Each offset measurement has an 
associat ed un certainty (a/?) dependent on atie, ctag, and ahost 
(Section [3T] ). Due to the inherent non-Gaussian nature of off¬ 
sets and the non-uniform uncertainties, we employ a Monte 
Carlo approach to assess the uncertainty on the resulting dis¬ 
tributions of offsets. At each iteration a random offset (x) was 
drawn for each burst from its offset probability distribution 
deflned by the measured offset (R) and uncertainty (a/?). For 
bursts with R/ag < 5 we use a Rice distribution to represent 
the probability distributions, deflned by: 


p{x\R,aR) = —ycxp 


(x^+R^)] fxR\ 

J V<^i/ 


( 2 ) 


where Iq is the zeroth order modifled Bessel function of the 
first kind. For R/aR > 5, the Rice distribution can be approxi¬ 
mated as a Gaussian. We make this approximation for numer¬ 
ical ease. In all cases discussed in the following sections we 
use 10,000 iterations. 


4.1.1. Physical Offset Distribution 

In Figure we show the cumulative distribution of pro¬ 
jected physical offsets, calculated as described in Section 


3.1 The final sample size included in this distribution is 71 
LGRBs, the subset of the original sample where offset mea¬ 
surement is possible. The distribution ranges from offsets of 
0.075 to 14.4 kpc with the median at 1.42 kpc. About 90% of 
bursts occur within ^ 5 kpc. The distribution is remarkably 
smooth and has bursts representing the full range of offsets. 
There are few gaps where we do not And bursts, with one 
notable gap at ^ 7 -14 kpc. For comparison we also show 
the distribution of ph ysical offsets for the small sample of 20 
LGRBs measured by [Bloom et al. ( |2002| ). The two distribu¬ 
tions agree exceptionally well. A Kolmogorov-Smirnov (KS) 
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Figure 5. Left: Cumulative distributions of the projected physical offsets of 71 LGRB s from this work (red) andjBloom et al.U2002| (B02, blue). We also show 
the distributions for Type Ib/c (cyan) and Type II (green) SNe fromjPrieto et al.U2008f The shaded region displays the results of our Monte Carlo assessment of 
the uncertainties on the offsets in the form of a 2D histogram. Right: Similar to the left panel, except here we plot the sample of 58 bursts with ^1-0 kpc 
to reduce any bias caused by bursts with a large offset uncertainty (See Figurej^. The results are essentially unchanged. 


test yields a p-value of 0.94. We also show in Figurethe re¬ 
sults of our Monte Carlo simulation (taking into account the 
offset uncertainties) by plotting a 2D histogram showing the 
density of points from the resulting cumulative distributions 
generated at each iteration. Darker regions indicate a higher 
density of points, in other words a high fraction of synthetic 
distributions pass through that region. The apparent shift of 
the distribution from the measured distribution towards higher 
offset values is indicative of the skew to higher offsets result¬ 
ing from the fact that an offset is a positive-definite quantity 
described by a Rice distribution. The shift is dominated by 
the ^10 LGRBs for which the uncertainty is larger than the 
offset. 

It is important to carefully consider the effects of the un¬ 
certainties on the distribution of offsets because offsets with 
uncertainties large relative to the offset itself have probabil¬ 
ity distributions skewed to larger offset. Given a large uncer¬ 
tainty, the measured offset is more likely to be large due to the 
large area encompassed by the error circle at large offset. As 
should not be correlated with the intrinsic offset, if the 
true offset distribution extended to small offsets, this effect 
would cause a bias against finding them when the uncertain¬ 
ties are large. In Figure]^ we plot cr/?pj^y^ versus Rphys- Of the 
13 bursts in the sample with ^1-0 kpc, all are at mea¬ 
sured offsets ^0.5 kpc, whereas the remaining bursts with 
^^phys ^ kO kpc span the full range of offsets from 0.075 to 14 
kpc. Since there is no relationship between the true offset and 
, the range of offsets found with large uncertainty should 
not be significantly different from the range found with small 
uncertainty. The fact that we only measure large offsets when 
we have large uncertainty is indicative of the above bias. We 
note that even at relatively small uncertainty we may still be 
missing bursts at very small offsets because any amount of 
uncertainty prevents the measurement of arbitrarily small off¬ 
sets. 


To reduce this bias, we make a quality cut on our physi¬ 
cal offset sample by also plotting in Figure the distribution 
for bursts with ^1-0 kpc, resulting in a sample of 58 
LGRBs. The resulting distribution has a median at 1.27 kpc. 
The mean of the distribution of medians for each iteration of 
the Monte Carlo simulation is 1.34 kpc with a 90% confidence 
interval of 1.17 - 1.51 kpc. 

We also compare the distribution of projected physical off¬ 
sets to the distributions for Type II and Ib/c SNe . KS tests with 
the supernovae samples of|Prieto et al.|(|2008|) yield p-values 
of 1.33 X10“® and 1.53 x ir'for the Type Ib/c and II SNe, re- 
spectively. LGRBs occur significantly closer to the centers of 
their host galaxies in a physical sense than SNe. To disentan¬ 
gle the effects of galaxy size differences between LGRB hosts 
and SNe hosts, we revisit supernova comparisons in the next 
section where we analyze the LGRB host-normalized offsets. 

4.1.2. Host-Normalized Offset Distribution 

Normalizing the offsets by the host galaxy sizes enables a 
fair comparison across the LGRB sample and with the SNe 
samples. In Figure |7] we plot the cumulative distribution of 
host-normalized offsets {R/Rk). The distribution ranges from 
about 0.04 to 7.0. The median of the distribution is 0.79 and 
about 90% of the bursts occur within a host-normalized offset 
of 2.2. The distribution is overall reminiscent of an expo¬ 
nential disk profile, the expected surface brightness profile of 
star forming disk galaxies. For comparison we also plot the 
predicted distribution of host-normalized offsets if LGRB lo¬ 
cations exactly trace an exponential disk profile. Although 
the shapes are overall consistent between our measured dis¬ 
tribution and that of an exponential disk, there is a notable 
shift to lower offsets seen in our distribution. Again, to com¬ 
pare with the previous offset sample of LGRBs w e show the 
host-normalized offsets from |Bloom et al.| ( |2002| ). A KS test 
between this sample and our sample gives a p-value of 0.50. 
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Figure 6. Left: Uncertainty in physical offset versus physical offset (/^phys)- The green line marks = /^phys- Right: Uncertainty in host-normalized 

offset {cfr^q^) versus host-normalized offset (Rnormf Similarly, the green line marks cfr^^^ = /^norm- Since cfr^^^ and are not related to the intrinsic offset, 
the apparent lack of small offsets with large uncertainty (the gap in the upper left quadrant of both plots) is indicative of a bias to large offset associated with 
bursts that have poor localizations. The dashed blue lines at = 1 kpc and = 0.5 mark the locations of our cuts. 


As a result of the small sample size, |B loom et al.|(|2002|) con¬ 
cluded that the LGRB host-normalized offsets were consistent 
with being drawn from the exponential disk profile. With our 
larger sample we can rule this out at a significance level of 
0.01. Given the irregular morphologies of LGRB host galax¬ 
ies this result is not surprising even if LGRBs trace star for¬ 
mation. 

As with the physical offsets, we also show in Figure [7] the 
uncertainty region determined from our Monte Carlo proce¬ 
dure. Again, we carefully consider the effects of the uncer¬ 
tainties on the distribution of host-normalized offsets to avoid 
the aforementioned bias associated with the offset measure¬ 
ment of bursts with large offset uncertainties. In Figure we 
also plot versus Rnorm- Of the 20 bursts in the sample 
with (Jr ^,0.5, all are at measured offsets > 0.4, whereas 
the remaining bursts with ^0-5 span the full range of 
offsets from 0.04 to 7.0. As with the physical offsets, our 
inability to measure small offsets when we have large uncer¬ 
tainty is indicative of a bias. Again, it is likely the case that 
even at relatively small uncertainty we may still be missing 
bursts at very small offsets. 

As before, we make a quality cut on our host-normalized 
offset sample by restricting the sample to bursts with (Jr^^^^ < 
0.5, resulting in a sample of 52 LGRBs. The resulting dis¬ 
tribution has a median at 0.54. In Figure we show the in¬ 
dividual probability distributions of each LGRB used in the 
Monte Carlo simulation. The mean of the distribution of me¬ 
dians for each iteration of the Monte Carlo simulation is 0.67 
with a 90% confidence interval of 0.58 - 0.76. We plot the 
distribution and Monte Carlo results in Figure A KS test 
with the exponential disk model yields a p-value of 2.6 x 10“"^. 
Furthermore, only ^ 4% of the Monte Carlo synthetic distri¬ 
butions have p-values >0.05 when compared to an exponen¬ 


tial disk profile, meaning that we can rule out the hypothesis 
that LGRB host-normalized offsets are drawn from an expo¬ 
nential disk. As another way to view these results we also 
show in Figure the sum of the individual probability distri¬ 
butions of each LGRB. The summed distribution effectively 
produces a smoothed histogram of the offsets and shows a 
clear shift to smaller offsets from the probability distribution 
for an exponential disk profile. LGRBs are apparently highly 
concentrated in the inner parts of their hosts. In particular, 
they are more concentrated than the radial light distributions 
of their hosts; namely, 50% of LGRBs occur within 67% of a 
half-light radius. Stated differently, about half of all LGRBs 
occur within a region that contains only 33% of the underly¬ 
ing distribution of light in their host galaxies. The possible 
bias caused by the inability to find bursts at very small offset 
mentioned above would act in the opposite direction of the 
striking result found here, indicating that the distribution may 
be even more concentrated than we find here. 

A striking feature of the distribution in Figure is the ap¬ 
parent tail to large offset extending from a gap at Rnorm ^ 
1.6-3.3. While it is not unexpected, given the rough simi¬ 
larity of our distribution to an exponential disk, to find bursts 
at large offset, it is somewhat surprising that ^10% of LGRBs 
appear to be located at Rnorm ^ 3. These bursts include GRBs 
050820, 081008, 080928, 060418, and 060505. We revisit the 
properties of these bursts and their assigned hosts. The host 
of GRB 050820 has a morphology consisting of a bright core 
with a diffuse tail with the GRB located on this diffuse emis¬ 
sion, leading to a host-normalized offset of 3.8 from the bright 
core. The burst’s location on the underlying galaxy emission 
and the Fee of 0-01 do not give reason to doubt the host as¬ 
sociation, the large offset being due to the unusual morphol¬ 


ogy. As discussed in Section |3.3.1 we assigned as the host 
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Figure 7. Cumulative distributions of LGRB host-normalized offsets for our sample (red) and the sample of |Bloom et al.M2002) (B02, blu e) shown on linear 
(left) and log (right) x-axis scales. Also shown are the distributions for Type Ib/c (cyan) and Type II (green) SJNe from|Kelly & KirshnerU2012| and SLSNe 
(orange) from|Lunnan et al.jpOlSf We also plot the distribution expected for an exponential disk profile (magenta). The shaded regions depict the uncertainties 
on our offset measurements determined using a Monte Carlo simulation. We find an offset distribution qualitatively similar to an exponential disk profile but 
statistically inconsistent. 




Figure 8. Left: Probability distributions for each LGRB with < 0.5 given its measured offset and offset uncertainty calculated using Equationl^ Narrow 

peaks indicate well-localized bursts while broader distributions indicate poorer localizations. Right: Sum (red) of the individual probability distributions (left) 
effectively producing a smoothed histogram of the offsets. In magenta we show the probability distribution for an exponential disk. The summed distribution is 
clearly offset to lower offsets than what would be expected for a population tracing an exponential disk. 


of GRB 080928 a non-coincident galaxy with a luminosity 
consistent with the observed distribution of LGRB host lu¬ 
minosities. However, we cannot rule out the possibility that 
the host was not detected and that the assigned candidate is in¬ 


stead the source of an intervening absor 


)tion system observed 
casting doubt on the 


_ ig absorp t] 

in the afterglow ( Vreeswijk et al.|2008| ), ( 
large inferred offset. GRB 081008, another burst with a host 
assigned based on our luminosity analysis, has a large normal- 
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Figure 9. Similar to Figure ^except here we show the distribution of host-normalized offsets only for bursts with < 0.5. The distribution is notably 

shifted to lower offset with a median of 0.54. This is because we have removed bursts with large offset uncertainties which bias the distribution to larger offsets 
(See Figurej^. 


ized offset of and Pec ~ 0.07. While the host assignment 
matches our criteria, we acknowledge the possibility that the 
assignment is not correct. For GRB 060418 we assign as the 
ho st the brighter of the t wo nearby host components identified 
by|Pollack et al.|(|20091). Although these components are dis- 
tinct from each other (see Figure^, the assignment of one or 
the other has little effect on the offset and no effect on the frac¬ 
tional flux. GRB 060505 is a controversial z = 0.089 burst that 
did not have an associated supernova as expected for a low-z 
LGRB leading som e to argue in favor of a short GRB scenario 
( [Qfek et ^|2007| ). However, given its location in a bright, 
low-metallicity HII region of its host, o thers argue in favor of 
060505 as a LGRB ( |Thone et al.|2014| ). In this case, 060505 
is an example of a LGRB occurring at large normalized off¬ 
set (3.6) and at high fractional flux (0.99) in the outskirts of 
its host. It is interesting to speculate that the other bursts at 
large offset may be higher redshift versions of 060505 where 
underlying HII regions are too faint to be detected. For the 
reasons discussed here, the large offset tail in the distribution 
of host-normalized offsets may not be real. However, it does 
not affect the overall conclusion that LGRBs as a population 
are more centrally concentrated than the underlying light dis¬ 
tributions of their hosts. 

We compare our distribution of LGRB ho st-normalized off¬ 
sets wi th those measured for supernovae by [Kelly & Kirshner 


( |2012| ), also shown in Figure KS tests with the Type Ib/c 
and Type II SNe yield p-values of 1.3 x 10“^ and 5.2 x 10“^, 
respectively. While the measured distribution for Type Ib/c 
and our sample of LGRBs are statistically inconsistent, about 
80% of Monte Carlo synthetic distributions have p-values 
> 0.05 when compared to Type Ib/c SNe. This percentage 
is ^1% when compared to Type II SNe. When taking into ac¬ 
count the uncertainties, the distributions for LGRBs and Type 
Ib/c SNe are therefore consistent. They both are much more 
centrally concentrated than Type II SNe and the exponential 
disk profile. Comparing our LGRBs to the distribution mea- 


sured for super-luminous supernovae (SLSNe; [Lunnan et al. 
2015| ) we find a KS test p-value of 0.17. This comparison, 
though, is limited by the small SLSNe sample size. 

4.1.3. Offset Distribution Summary 

We find that our sample of LGRBs spans the range 0.075 
to 14.4 kpc (mean of Monte Carlo distribution of medians 
= 1.34 kpc) in physical offset from the centers of their host 
galaxies. When appropriately accounting for uncertainties 
and removing bursts with ^0-5, we find that LGRB 
host-normalized offsets are considerably smaller and more 
centrally concentrated (mean of Monte Carlo distribution of 
medians = 0.67, with 90% confidence interval of 0.58 - 0.76) 
than what would be expected if LGRBs traced an exponential 
disk profile. We can rule out at high significance that LGRB 
offsets are drawn from an exponential disk distribution and 
their distribution is inconsistent with Type II SNe. Account¬ 
ing for uncertainties, the LGRB offset distribution we find is 
consistent with that for Type Ib/c SNe. LGRBs prefer the cen¬ 
tral locations in their hosts, being more concentrated than the 
underlying host light distributions. Table provides a sum¬ 
mary of the KS p-values we find for comparisons between 
our sample of LGRB host-normalized offsets and other types 
of stellar explosions and in Tablewe show the medians. 

4.2. Host Light Distribution 

In Figure we show the cumulative distribution of frac¬ 
tional fiux for the 49 bursts in our sample with error circle to 
galaxy area ratios of <0.1 (left panel). Ideally, we would 
like to assess whether or not LGRBs are spatially coinci¬ 
dent with bright rest-frame UV regions of their hosts, which 
trace young, massive stars. As shown in Figure our ob¬ 
servations primarily probe the rest-frame optical emission, 
which may probe older stellar populations. Here we com¬ 
pare the fractional fiux distribution for the ^ 20% of the sam¬ 
ple that probes rest-frame UV to the distribution that probes 
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rest-frame optical. The right panel of Figureshows the re¬ 
sulting distributions when dividing the sampleoy wavelength 
into two bins separated at Aj-est = 5000 A. The rest-frame opti¬ 
cal and blue/UV distributions agree well, which we interpret 
as indicating that there is no bias in using rest-frame opti¬ 
cal observations as a proxy for rest-frame UV. If there was a 
separate older stellar population, we would not expect to find 
LGRB positions correlated so strongly with both the optical 
and UV light. In other words, the optical and UV light are 
probing the same population of stars. We can then interpret 
our fractional flux measurements as refiecting the relationship 
between location and the host UV emission. 

As shown in Figure the fractional flux distribution for 
our sample spans the fuITrange from 0 to 1 with 90% of val¬ 
ues falling in the range 0.03 to 0.98. The median is 0.75 and 
nearly 80% of bursts have fractional flux values greater than 
0.5. This strong preference towards higher fractional flux val¬ 
ues indicates that LGRBs are typically located on some of the 
brightest regions of their hosts, and lends further support to 
the idea that LGRBs originate from very mas sive stars. This 
is in broad agreement with previous claims ( [Fruchter et al. 
2006} [Svensson et al.||201Q| ). We also show in Figure |10| ^ 
distribution of fractional fiux for the sample of 30 LGfi 


tEe 

RBs 


_ sample 

from [Fruchter et al.] ( |2006| ). A KS test between the two sam¬ 
ples yields a p- value of 0.09. This i ndicates that while our 
sample and the Fruchter et al. ( 2006| ) sample are marginally 
inconsistent, we cannot rule out that they are drawn from the 
same d istribution. Nevertheless, the sample of [Fruchter et al.j 
( [2006| ) has a higher median and general inclination towards 
higher fractional fiux values. This slight discrepancy is prob¬ 
ably not caused by effects intrinsic to the LGRB population 
but rather by syst ematic differences be tween the two samples. 
We also show the [Fruchter et al.[ (2006 ) sample with the wave¬ 
length binned distributions in Figure 1^ The slight shift of 
th e rest-frame blue/UV distribution towards the distribution 
of [Fruchter et al.[p006[ ), mostly comprised of rest-frame UV 
observations, may indicate the discrepancy is a bandpass ef¬ 
fect. We address another potential effect in Section 

We also compare our measured distributi on to the corre- 
sponding distributions measured for CCSNe ([Svensson et al. 
[201()[), SLSNe (ILunnan et al.[[2015| ), and Type Ic and II SNe 
( [Kelly et al.|2008[ ). KS tests between our sample, the SLSNe 
sample, the Type Ic SNe sample, the Type II SNe sample, and 
the CCSNe sample yield p-values of 0.52, 0.79, 5.4x10“^, 
and 0.03, respectively. We cannot rule out that our sample of 
LGRBs and the SLSNe sample are drawn from the same dis¬ 
tribution. Interestingly, our sample overlaps the SLSNe distri¬ 
bution at fractional fiux values of > 0.8. The very good agree¬ 
ment between the distributions for our sample and the Ic SNe 
sample is consistent with the fact that Ic SNe, particularly Ic- 
BL, are the only types of SNe observed to be associated with 
LGRBs. The distributions of fractional fiux for the CCSNe, 
mostly consisting of Type II SNe, and the Type II SNe sam¬ 
ple are clearly distinct from our LGRB distribution, strong 
evidence that Type II SNe and LGRBs have different progen¬ 
itors. We summarize the KS p-values for our comparisons in 
TablelH 

4.3. Galaxy Sizes 

As it pertains to our knowledge of the host galaxies and 
therefore implications for the progenitors, we plot in Fig¬ 
ure [^ the cumulative distribution of 80% light-radii (Rgo) 


pie is 3.0 kpc and the distribution spans from Rgo ^ 0-6 kpc to 
^80 10 kpc. For comparison, we show the distributions for 

the sam ples of LGRBs and CCSNe studied by [Svensson et al.[ 
( 2010| ). A KS test reveals that our distribution, though shifted 
to slightly higher values of Rgo. is consistent with their dis¬ 
tribution of LGRB sizes. A KS test between our distribution 
of LGRB sizes and the CCSNe distribution yields a p-value 
of 3 X 10“^ indicating our larger sample of LGRBs is in full 


agreement with the conclusion of [Svensson et al.[ ( |2010[ ) that 
LGRB host galaxies are on average smaller than CCSNe host 


galaxies. We also show in Figure 

the sample of SLSNe 

host galaxies studied by[Lunnan et al. 

|(2015[). With the pre- 

vious small sample of LGRBs, the SLSNe and LGRB host 


of our host galaxies. The median size for the LGRB sam- 


size distributions were statistically consistent. Now with our 
larger sample we find that the apparent larger sizes of LGRB 
hosts is statistically robust (KS p-value = 1.4 x 10“^). 

4.4. Durations and Energies 

The durations (as measured by T^q) and isotropic equivalent 
energies (£’ 7 ,^ 0 ) of LGRBs each vary by over several orders 
of magnitude. Whether or not environmental factors influence 
these variations in burst properties is an open question that we 
can address with our sample of offsets and fractional fluxes. 
In Figure we plot fractional fiux and host-normalized off¬ 
set versus £90 and We find no significant correlations 

between fractional fiux or offset with £90 or E^^iso, indicating 
that the locations of LGRBs within their host galaxies do not 
affect the properties of the bursts themselves. Our findings 
do not support the claim of a possible trend between host- 
normalized offset and isotro pic equivalent energy reported by 
Ramirez-Ruiz et slI] ( [2QQ2| ) based on a much smaller sample 
of 16 prc-Swift events. 

5. DISCUSSION 

The distributions of offsets and fractional fiux presented in 
Section [^provide the most in-depth look at the locations and 
environments of LGRBs presented to date. Furthermore, our 
sample is sufficiently large and wide in redshift to allow inves¬ 
tigation of trends with cosmic time. Here we discuss possible 
redshift evolution, the relationship between fractional fiux and 
host-normalized offset, and the implications for the progeni¬ 
tors of LGRBs provided by this new view of their locations. 

5.1. Redshift Trends 

In Figure 13 we show the physical offsets for bursts with 
^^phys ^ LO^c as a function of redshift (right panel) and 
the cumulative distribution of physical offsets binned into 
three equally populated bins of redshift (left panel): z < 0.9, 
0.9 < z < 1. 8 , and z > L 8 . We find that there is no statistically 
significant trend of the physical offsets as a function of red¬ 
shift. KS tests between the three redshift bins yield p-values 
of 0.97, 0.77, and 0.45 for the low-z and mid-z bins, low-z and 
high-z bins, and mid-z and high-z bins, respectively. 

In Figure 14 we show the host-normalized offsets of bursts 
with <TJ;5as a function of redshift (right panel) and the 
cumulative distribution of host-normalized offsets binned by 
redshift (left panel). Overall we find no significant trend of 
the host-normalized offsets with redshift. KS tests between 
the three redshift bins yield p-values of 0.61, 0.51, and 0.23 
for the low-z and mid-z bins, low-z and high-z bins, and mid-z 
and high-z bins, respectively. Despite the statistical insignif¬ 
icance, there appears to be a slight shift of the z > L 8 bin 
toward larger offsets. If this reflects an actual trend, coupled 
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Figure 10. L eft: Cumulative dist ributions of fractional flux for our sample (red), the LGRB sample of|Fruchte r et al.H2006) (F06 , green), CCSNe from|Svensson| 
|et al.|j2010) (blue), SLSNe from|Lunnan et al.|J2015| (orange), and Type Ic (cyan) and II (violet) SNe trom |Kelly et ^ 2008). We only include bursts m our 
sample with a ratio of error circle to galaxy area of ^ 0.1 (Figure]^. The diagonal dashed line shows the expected tractional flux distribution for a population 
uniformly tracing the underlying light of their host galaxies. Right: tumulative distributions of fractional flux for two subsets of the sample divided at a rest-frame 
wavelength of 5000 A. The agreement between the two means that we can use the rest-frame optical observations as a proxy for rest-frame UV. 



Figure 11. Cumulative distributions of Rgo for our s ample of LGRBs (red), 
the LGRB (green) and CCSNe (bl ue) samples fromlSvensson et al.H2010) 
(SIO), and the SLSNe sample from|Lunnan et al.|fl2015) (orange). At high 
signiflcance we And that LGRB host galaxy sizes are larger than SLSNe host 
sizes but smaller than CCSNe host sizes. 


z < 2 the median galaxy half-light radius is 1.98 kpc and at 
z > 2 the median radius is 1.15 kpc. It is possible that sur¬ 
face brightness dimming, a strongly redshift dependent effect 
((X (1 +z)“^), may be limiting our ability to measure the sizes 
of the host galaxies at higher redshift. At higher redshifts the 
faint outer regions of galaxies become much harder to detect 
resulting in an underestimate of galaxy size and correspond¬ 
ing overestimate of the host normalized offset. However, HST 
observations of Lyman-break galaxies indicate that high red¬ 
shift galaxies are intrinsically co mpact and show a size evo¬ 
lution of ap proximately (1 +z)~^ ( Bouwens et al.|2004 


lution ot app 
et al.|j2014| ). 

lack diffuse ( 


Calvi 


If low-mass galaxies such as LGRB hosts also 
outer components at high redshift, then it is possi¬ 
ble the observed size trend in Figure [T^ reflects actual galaxy 
growth. 

In Figure[^we show the fractional flux distribution binned 
into the same three redshift bins. KS tests between the three 
redshift bins yield p-values of 0.92, 0.13, and 0.42 for the 
low-z and mid-z bins, low-z and high-z bins, and mid-z and 
high-z bins, respectively. While the three binned distributions 
are statistically consistent, the lowest fractional flux values (< 
0.2) only occur at z > 1.5. It is possible that surface brightness 
dimming is also playing a role here. 

Overall we find no statistically significant trends with red¬ 
shift in the offset or fractional flux distributions. The slight 
shift towards larger offsets and lower fractional flux in the 
highest redshift bin of z > 1.8 may indicate we are missing 
some of the diffuse outer parts of the hosts at high redshift. 


with the lack of a trend in the distribution of physical offsets, 
it would indicate an effect related to the sizes of the galax¬ 
ies. In Figure we show the measured half-light radii of 
the host galaxies as a function of redshift. There is clearly 
a trend of smaller galaxy sizes with increasing redshift. At 


5.2. The Fractional Flux - Offset Relationship 

In Figure we plot the fractional flux versus host- 
normalized oteet for bursts satisfying error circle area to 
galaxy area ratio of <0.1 and ^0.5. We find a clear 
correlation between fractional flux and host-normalized off- 
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Figure 12. Scatter plots of fractional flux and host-normalized offset {Rnomd versus duration (T^q) and equivalent isotropic energy (E^^iso)- We And no significant 
correlations between these quantities. 




Figure 13. Left: Cumulative distribution of physical offsets for bursts with ^1-0 kpc binned into three equal-numbered bins of redshift: z < 0.9 (blue), 
0.9 < z < 1.8 (red), and z > 1.8 (green). Right: Physical offsets with their associated uncertainties as a function of redshift. In both plots we do not include bursts 
with unknown redshift. We do not find a statistically significant trend in physical offset as a function of redshift. 
































20 


Blanchard et al. 





a; 

i/) 

tt 

O 

■o 

(D 

_N 

75 

E 


{/) 

o 

X 



2; 


Figure 14. Same as Figure [T3]but for the host-normalized offsets of bursts with < 0.5. We do not find a statistically significant trend in host-normalized 

offset as a function of redshim The slight inclination towards larger offsets at higher redshift may be the result of surface brightness dimming. 



Figure 15. Half-light radii as a function of redshift showing a clear trend 
that may explain the slight shift to larger host-normalized offsets at z > 2. 
Whether this trend reflects real galaxy growth or the effects of surface bright¬ 
ness dimming is less clear. 


set, where bursts at smaller offsets have high fractional flux 
values and bursts at larger offsets have on average lower frac¬ 
tional flux values. The scatter in the fractional flux values 
increases with increasing /?norm- All bursts in our sample with 
^norm ^0-5 have fractional flux values greater than 0.6, yield¬ 
ing a notable lack of bursts with small offsets and low frac¬ 


tional flux values. Beyond /?norm ~ 0.5 the fractional flux val¬ 
ues are on average lower but also show considerable scatter. 
Bursts with Rnonn ^ 1 have fractional flux values spanning the 
full range from 0 -1. 

While it is not surprising that a trend should exist between 
fractional flux and offset given that regions close to the centers 
of galaxies are brighter, the coupling of fractional flux and off¬ 
set information explains the results seen in the fractional flux 
distribution. To illustrate this, we also show in Figure the 
cumulative distributions of fractional flux for the subsamples 
with 7?norm < 0-5 and Rnovm > 0.5. Bursts with /?norm < 0.5 
exclusively occur on regions of high fractional flux and bursts 
with Rnovm >0.5 uniformly trace the light of their hosts, show¬ 
ing no preference for unusually bright regions. This means 
that the shift to high fractional flux we And in the total dis¬ 
tribution is entirely due to LGRBs at small offsets and not to 
bursts at large offset that happen to reside in bright star form¬ 
ing regions. Therefore, it is not simply the presence of a bright 
star forming region that affects LGRB production, but rather 
another factor playing a role near the centers of galaxies must 
be at play. 

Previous studies indicated that the agreement between 
LGRB offsets and the exponential disk was inconsistent with 
the observation that LGRBs we re highly correlated with the 
brightest re gions of their hosts ( [Bloom et al.||2QQ2t [Fruchterj 
et al.|2006] ). The offset distribution we measure coupled with 


the cumulative distribution of fractional flux showing a pref¬ 
erence for bright regions does not imply a discrepancy. In¬ 
deed, the fact that we And LGRBs concentrated in the inner 
regions of their galaxies with corresponding high fractional 
fluxes can be i nterpreted as a relief of t he tension between th e 
conclusions of [Bloom et al.| ( |20Q2| ) and [Fruchter et al.[(2006[ ). 
Ultimately, the slight difference betw een our distribution of 
fractional flux and the distribution of [Fruchter et al.[ ( [2006[ ) 
may be explained by a difference in the offset distributions of 




















































An HST Study of the Locations of Long Gamma-Ray Bursts 


21 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 


Figure 16. Same as FigurefT^but for fractional flux. We do not And a statistically signiflcant trend in fractional flux as a function of redshift. 




Figure 17. Left: Fractional flux vers us h ost-normalized offset for L GRBs (red) with ^ 0.5 and error circle to galaxy area ratio <0.1 and Type Ic/Ic-BL 

SNe (cyan) from [Kelly et al.U2008) and [Kelly & Kirshner|(|2012). Right: Cumulative distributions of fractional flux for two subsets of the LGRB sample 
divided at Rnorm = 0.5. Ihere is a clear trend relating these two quantities where LGRBs at small host-normalized offset, near the bright cores of their hosts, are 
exclusively at high fractional flux. At Rnorm ^ 0.5, the bursts uniformly trace the underlying light of their host galaxies similar to CCSNe as also shown in the 
right panel. Therefore, the shift to high fractional flux seen in the full sample distribution (Figure is due to bursts occurring in the inner parts of their host 
galaxies. A similar trend is seen in the Ic SNe sample. 


the two samples, and may simply be the result of small num¬ 
ber statistics where our sample more closely approximates the 
true underlying distribution. 


In Figure we also plot the normalized offsets and frac¬ 
tional fluxes m Type Ic and Ic-BL SNe for objects common to 
both the sample with fractional fluxes from|Kelly et al.|(|2008 ) 
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and the sample with host-normalized offsets from [Kelly & 
[Kirshner] ( |2012| ). We find a similar trend between the two 
quantities that indicates the shift to high fractional fiux in the 
Ic SNe distribution can be attributed to SNe that occurred at 
small offset. 


5.3. Progenitor Implications 

We find that LGRBs are on average located at smaller host- 
normalized offset than expected for an exponential disk pro¬ 
file of star formation and are correspondingly more concen¬ 
trated on brighter regions of their hosts than expected for a 
population simply tracing the underlying light distributions. 
The evidence presented here suggests that the progenitors of 
LGRBs prefer the central regions of their hosts, indicating 
that the properties of star formation towards the inner regions 
of LGRB host galaxies are favorable for LGRB production. 
As demonstrated in Figure [T7] the bursts responsible for the 
strong preference towards mgher fractional flux values are 
dominated by bursts at small offsets. The fractional flux dis¬ 
tribution of bursts at offsets Rnorm ^0-5 uniformly tracks the 
host light. In effect the high fractional fiux values are due to 
the fact that the bursts occur at small offsets where the host is 
brighter. 

How the suspected metallicity bias of LGRB production fits 
into this new view of their preferred locations is an important 
question. Star forming disk galaxies typically have me tallic 
ity gradients with higher metallicities towards th e cores ( [Vila 
[Costas & Edmunds] 1992[ [Zaritsky et al.[[1994| ). One might 
then expect, if metallicity gradients are also present in the 
hosts of LGRBs, to more frequently find LGRBs in the out¬ 
skirts of their hosts contrary to what we find here. Recent re¬ 
solved host stud ies of low redshift hosts indicate weak m etal- 
licity gradients ( [Levesque et aH20 fT] [Thone et al.|2014| ) and 
that the global host metallicity is often a sufficient proxy for 
the metallicity of the burst site. Though the metallicity typ¬ 
ically shows only slight variation a cross the host, the burst 
site often has the lowest metallicity ([Christensen et al.|2008[ 


[Levesque et al.|2QTT| [Thone et al.|26l4[ ). In future work, we 

plan to test for the presence of metallicity gradients, in an av¬ 
erage sense, in high redshift hosts by combing ISM metallic¬ 
ity measurements with our measured offsets. While it is clear 
that metallicity is a factor affecting the production of LGRBs, 
our results indicate that other factors relating to the mode of 
star formation may also play a dominant role. The conditions 
in the central regions of LGRB hosts, perhaps a consequence 
of their formation histories, may be such that the IMF is dif¬ 
ferent, producing more massive stars and thus more potential 
LGRB progenitors. It has been shown that massive stars are 
oft en fou nd in binary systems ( [Mason et al.[[2QQ9'} [Sana et[ 
[ar ||2012| ) and that binary interactions may be important fo r 
LGRB production ( [Fryer et al.|1999[[Bromm & Loeb|2006| . 
There may be increased numbers of massive star binaries in 
the central regions of LGRB hosts. 

In both the offset and fractional flux distributions we find 
that LGRBs and Type Ic SNe prefer similar locations, in¬ 
dicating a common link between their progenitors. While 
LGRBs and Type Ic SNe broadly prefer similar locations, the 
exclusive association of LGRBs with Ic-BL SNe indicates that 
some factor must determine whether a star explodes as a nor¬ 
mal Type Ic, a Type Ic-BL without a LGRB, or a Type Ic-BL 
with a LGRB. Previous observations suggest that the environ¬ 
ments where Type Ic-BL SNe occur are more metal poor than 
those for normal Type Ic SNe, and that the environments of 
Type Ic-BL SNe associated with LGRBs are even more metal 


poor ( [Modjaz et al.|2008l[20111 Kelly & Kirshner|2012| ). Re 


cent work suggests that nikallicity is likely not the only factor 
infiuencing the production of Ic-BL SNe and LGRBs. [Kelly[ 
[et al.[ ( |MT4l ) found that Ic-BL SNe and LGRBs occur in host 
galaxies that have higher stellar mass and star formation rate 
densities than a control sample of SDSS star forming galax¬ 
ies. They show that metallicity cannot explain this preference 
and instead suggest that another environmental factor such as 
increased numbers of young, bound clusters where tight bina¬ 
ries form or a top-heavy IMF may be at play in Ic-BL SNe and 
LGRB hosts. The results we find, namely that LGRBs prefer 
the central regions of their hosts, may be a sub-galactic mani¬ 
festation of the same environmental factor responsible for the 
preference for hi gh global stellar m ass and star formation rate 
densities seen by Kelly et aL[ ( [2014[ ). Whatever this factor is, it 
seems to infiuence primarily where LGRBs occur but not the 
properties of the explosions themselves as we do not find cor¬ 
relations between fractional flux or offset with burst duration 
or isotropic equivalent energy. 

Although the SLSN sample size is much smaller than the 
LGRB sample size, the statistical consistency between their 
distributions of fractional fiux may be hinting that SLSNe and 
LGRBs share a common mass progenitor. It could be that 
in the case of SLSNe, the core is not quite massive enough 
to form a black hole and instead collapses into a magnetar, 
while in the case of LGR Bs, the core collapses into a black 
hole ( [Lunnan et al.|2014| ). However, the difference in the size 
distributions of SLSN and LGRB hosts indicates other factors 
also infiuence the production of their progenitors. 


6. CONCLUSIONS 

We have carried out a comprehensive study of the loca¬ 
tions of LGRBs within their host galaxies using HST observa¬ 
tions of Swift bursts compiled over the last decade. Using as¬ 
trometry from ground- and space-based detections of optical, 
NIR, radio, and X-ray afterglows, we measure the projected 
physical and host-normalized offsets of 71 LGRBs from their 
host centers. In addition, we measure the fractional flux, the 
brightness of the burst site relative to the total host light distri¬ 
bution. Upon restricting the sample to avoid biases associated 
with measurement uncertainty we obtain sample sizes of 58, 
52, and 50 LGRBs for the physical offsets, host-normalized 
offsets, and fractional fluxes, respectively. These measure¬ 
ments have enabled us to study the offset and fractional flux 
distributions in great detail, providing an in-depth view into 
the preferred locations of LGRBs. Our results are as follows: 


In ag reement with previous work ( [Wainwright et al. 


2007), we find that LGRB hosts typically show irreg 


ular morphologies and are on average smaller (median 
= 3.0 kpc) than the host galaxies of CCSNe, though not 
as small as the hosts of SLSNe. 


• We find that the physical offsets of LGRBs have a mean 
of the Monte Carlo medians of 1.34 kpc with a 90% 
confidence interval of 1.17-1.51 kpc. 

• The host-normalized offsets span a wide range of 0.04 
to 7.0 but are concentrated at small values with a mean 
of the Monte Carlo medians of 0.67, with a 90% con¬ 
fidence interval of 0.58 - 0.76. This distribution of 
LGRB offsets is considerably more concentrated than 
expected if LGRBs traced an exponential disk (p-value 
= 2.6 X 10“"^). LGRBs are more concentrated than their 
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host galaxies’ own light distributions, such that 50% of 
bursts reside within a region containing only 33% of the 
underlying host light. 

• Compared to Type Ib/c SNe, LGRBs have consistent 
host-normalized offset distributions with 80% of syn¬ 
thetic Monte Carlo distributions yielding p-values > 
0.05. On the other hand, LGRB offsets are inconsis¬ 
tent with those of Type II SNe. 

• The fractional flux distribution of our sample of LGRBs 
has a median at 0.75 and ^ 80% of bursts have frac¬ 
tional flux values ^0.5, indicating a strong preference 
towards high fractional flux though not as strong as pre¬ 
vious studies of smaller samples indicated. 

• The host-normalized offsets and fractional fluxes are 
correlated, such that bursts at small offset (< 0.5) occur 
exclusively on regions of high fractional flux (> 0.6) 
and bursts at larger offset uniformly track the underly¬ 
ing host light distributions, indicating that the prefer¬ 
ence for high fractional flux present in the full sample 
is entirely due to LGRBs at small offset. 

• The offset and fractional flux distributions show no sta¬ 
tistically signiflcant trends with redshift, though a slight 
shift towards larger offsets and lower fractional flux val¬ 
ues at z > 1.8 may indicate the effects of surface bright¬ 
ness dimming. 

From our study it appears that the star formation occurring 
within ^ Rh is the most favorable for LGRB production. The 
distributions of offset and fractional flux for our sample of 
LGRBs agree well with those for Type Ic SNe, a scenario 
consistent with the well-known connection between LGRBs 
and Type Ic-BL SNe. Both types of transients occur in spe¬ 
cial types of galaxies and in special environments within their 
hosts, where metallicity is likely not the only important factor. 
In the central regions of galaxies it is possible that changes in 
the IMF or increased massive star binary fractions may be re¬ 
sponsible for the sub-galactic preferences of LGRBs and Ic- 
BL SNe. This new view of LGRB locations we have pre¬ 
sented is unlikely to change in the near future because the 
sample size is unlikely to increase signiflcantly. Progress to¬ 
ward understanding the precise factors at play in the central 
regions of LGRB host galaxies will likely come through stud¬ 
ies of LGRB progenitor analogs in the local universe whose 
environments can be studied in great detail. 

The Berger GRB group at Harvard is supported in part by 
the NSF under grant AST-1411763 and by NASA under grant 
NNX15AE50G. This paper is based upon work supported by 
the National Science Foundation Graduate Research Fellow¬ 
ship Program under Grant No. DGEl 144152. W.F. acknowl¬ 
edges support provided by NASA through Einstein Postdoc¬ 
toral Fellowship grant number PF4-150121 issued by the 
Chandra X-ray Observatory Center, which is operated by the 
Smithsonian Astrophysical Observatory for NASA under con¬ 
tract NAS8-03060. We thank Ragnhild Lunnan, Maria Drout, 
Josh Grindlay, and Avi Loeb for helpful discussions. The 
HST data presented here and the SwiftfCWOT data used in 
this work were obtained from the Mikulski Archive for Space 
Telescopes (MAST). STScI is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA 


contract NAS5-26555. Support for MAST for non-HST data 
is provided by the NASA Office of Space Science via grant 
NNX09AF08G and by other grants and contracts. Based in 
part on observations obtained at the Gemini Observatory, ac¬ 
quired through the Gemini Science Archive and processed 
using the Gemini IRAF package, which is operated by the 
Association of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United 
States), the National Research Council (Canada), CONICYT 
(Chile), the Australian Research Council (Australia), Min- 
isterio da Ciencia, Tecnologia e Inova 9 ao (Brazil) and Min- 
isterio de Ciencia, Tecnologia e Innovacion Productiva (Ar¬ 
gentina). Based in part on data obtained from the ESO Sci¬ 
ence Archive Eacility. This paper makes use of data gathered 
with the 6.5 meter Magellan Telescopes located at Las Cam- 
panas Observatory, Chile. This research has made use of data 
obtained from the Chandra Data Archive and software pro¬ 
vided by the Chandra X-ray Center (CXC) in the application 
package CIAO. Based in part on data collected at the Subaru 
Telescope and obtained from SMOKA, which is operated by 
the Astronomy Data Center, National Astronomical Observa¬ 
tory of Japan. We thank Brad Cenko for providing Palomar 
60inch afterglow images for GRBs 060206 and 090618. 

REFERENCES 


Alard, C. 2000, A&AS, 144, 363 

Beckwith, S. V. W., Stiavelli, M., Koekemoer, A. M., et al. 2006, AJ, 132, 
1729 

Berger, E. 2010, ApJ, 722, 1946 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Bloom, J. S., Kulkarni, S. R., & Djorgovski, S. G. 2002, AJ, 123, 1111 
Bouwens, R. J., Illingworth, G. D., Blakeslee, J. R, Broadhurst, T. J., & 
Franx, M. 2004, ApJ, 611, LI 

Bouwens, R. J., Illingworth, G. D., Oesch, R A., et al. 2015, ApJ, 803, 34 
Bromm, V., & Loeb, A. 2006, ApJ, 642, 382 

Calvi, V., Stiavelli, M., Bradley, L., Rizzella, A., & Kim, S. 2014, ApJ, 796, 
102 

Christensen, L., Hjorth, J., & Gorosabel, J. 2004, A&A, 425, 913 
Christensen, L., Vreeswijk, R M., Sollerman, J., et al. 2008, A&A, 490, 45 
Faber, S. M., Willmer, C. N. A., Wolf, C., et al. 2007, ApJ, 665, 265 
Frail, D. A., Kulkarni, S. R., Sari, R., et al. 2001, ApJ, 562, L55 
Fruchter, A. S., Levan, A. J., Strolger, L., et al. 2006, Nature, 441, 463 
Fryer, C. L., Woosley, S. E., & Hartmann, D. H. 1999, ApJ, 526, 152 
Gonzaga, S., & et al. 2012, The DrizzleRac Handbook, HST Data Handbook, 
Graham, J. E, & Fruchter, A. S. 2013, ApJ, 774, 119 
Hogg, D. W., Rahre, M. A., McCarthy, J. K., et al. 1997, MNRAS, 288, 404 
Kelly, R. L., Kirshner, R. R, & Rahre, M. 2008, ApJ, 687, 1201 
Kelly, R. L., & Kirshner, R. R. 2012, ApJ, 759, 107 

Kelly, R. L., Filippenko, A. V., Modjaz, M., & Kocevski, D. 2014, ApJ, 789, 
23 

Kruhler, T., Greiner, J., Schady, R, et al. 2011, A&A, 534, A108 
Le Floc’h, E., Due, R.-A., Mirabel, I. E, et al. 2003, A&A, 400, 499 
Levesque, E. M., Kewley, L. J., Graham, J. E, & Fruchter, A. S. 2010a, ApJ, 
712, L26 

Levesque, E. M., Kewley, L. J., Berger, E., & Zahid, H. J. 2010b, AJ, 140, 
1557 

Levesque, E. M., Berger, E., Soderberg, A. M., & Chornock, R. 2011, ApJ, 
739, 23 

Levesque, E. M. 2012, Death of Massive Stars: Supernovae and 
Gamma-Ray Bursts, 279, 167 

Lunnan, R., Chornock, R., Berger, E., et al. 2014, ApJ, 787, 138 
Lunnan, R., Chornock, R., Berger, E., et al. 2015, ApJ, 804, 90 
MacFadyen, A. L, & Woosley, S. E. 1999, ApJ, 524, 262 
Marchesini, D., van Dokkum, R, Quadri, R., et al. 2007, ApJ, 656, 42 
Mason, B. D., Hartkopf, W. I., Gies, D. R., Henry, T. J., & Helsel, J. W. 
2009, AJ, 137, 3358 

Metcalfe, N., Shanks, T., Weilbacher, R. M., et al. 2006, MNRAS, 370, 1257 
Modjaz, M., Kewley, L., Kirshner, R. R, et al. 2008, AJ, 135, 1136 
Modjaz, M., Kewley, L., Bloom, J. S., et al. 2011, ApJ, 731, L4 



24 


Blanchard et al. 


Ofek, E. O., Cenko, S. B., Gal-Yam, A., et al. 2007, ApJ, 662, 1129 
Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 
Perley, D. A., Levan, A. J., Tanvir, N. R., et al. 2013, ApJ, 778, 128 
Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014, A&A, 571, 
A16 

Podsiadlowski, R, Mazzali, P. A., Nomoto, K., Lazzati, D., & Cappellaro, E. 
2004, ApJ, 607, L17 

Pollack, L. K., Chen, H.-W., Prochaska, J. X., & Bloom, J. S. 2009, ApJ, 
701, 1605 

Prieto, J. L., Stanek, K. Z., & Beacom, J. E. 2008, ApJ, 673, 999 
Ramirez-Ruiz, E., Lazzati, D., & Blain, A. W. 2002, ApJ, 565, L9 
Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444 
Savaglio, S., Glazebrook, K., & Le Borgne, D. 2009, ApJ, 691, 182 
Schlafly, E. E, & Finkbeiner, D. P. 2011, ApJ, 737, 103 
Stanek, K. Z., Gnedin, O. Y, Beacom, J. F., et al. 2006, AcA, 56, 333 


Svensson, K. M., Levan, A. J., Tanvir, N. R., Fruchter, A. S., & Strolger, 
L.-G. 2010, MNRAS, 405, 57 

Thone, C. C., Christensen, L., Prochaska, J. X., et al. 2014, MNRAS, 441, 
2034 

van den Heuvel, E. P. J., & Yoon, S.-C. 2007, Ap&SS, 311, 177 
Vila-Costas, M. B., & Edmunds, M. G. 1992, MNRAS, 259, 121 
Vreeswijk, P, Malesani, D., Fynbo, J., et al. 2008, GRB Coordinates 
Network, 8301, 1 

Wainwright, C., Berger, E., & Penprase, B. E. 2007, ApJ, 657, 367 

Woosley, S. E. 1993, ApJ, 405, 273 

Woosley, S. E., & Bloom, J. S. 2006, ARA&A, 44, 507 

Zaritsky, D., Kennicutt, R. C., Jr., & Huchra, J. P. 1994, ApJ, 420, 87 



An HST Study of the Locations of Long Gamma-Ray Bursts 

Table 1 

HST and Afterglow Observations 


HST Image Afterglow Image 


GRB 

z 

Instrument 

Filter 

Exp. Time 
(s) 

Obs. Date 
(UT) 

Program ID 

Telescope 

Filter 

040701 


ACSAVFCl 

F775W 

1940 

2004 Aug 09 

10135 

Chandra 





F625W 

1820 





040812 


ACSAVFCl 

F775W 

2120 

2004 Oct 04 

10135 

Chandra 





F625W 

2000 





040916 


ACSAVFCl 

F775W 

2056 

2004 Nov 30 

10135 






F625W 

1928 





040924 

0.859 

ACSAVFCl 

F775W 

3932 

2005 Feb 18 

10135 

/7^r/ACSAVFCl 

F775W 




F850LP 

3932 

2005 Feb 19 




041006 

0.716 

ACSAVFCl 

F775W 

4224 

2005 Feb 10 

10135 

/7^r/ACSAVFCl 

F775W 




F850LP 

4224 

2005 Feb 11 




050315 

1.949 

WFC3/IR 

F160W 

1209 

2011Jul 20 

12307 

Magellan/LDSS3 

r 

050401 

2.9 

WFC3/IR 

F160W 

1612 

2010 Oct 1 

12307 

VLT/FORS2 

R 

050406 

2.44 

WFC3/IR 

F160W 

1612 

2011 Feb 23 

12307 

Magellan/LDSS3 

r 

050408 

1.236 

WFC3/IR 

FI low 

1012 

2013 Mar 24 

12949 

Gemini-N/GMOS 

i 




F160W 

1209 





050416A 

0.654 

ACSAVFCl 

F775W 

4224 

2005 Nov 21 

10135 

7/*Sr/ACSAVFCl 

F775W 




F850LP 






050525 

0.606 

ACSAVFCl 

F625W 

4268 

2006 Mar 10 

10135 

//*Sr/ACSAVFCl 

F625W 




F775W 






050730 

3.967 

ACSAVFCl 

F775W 

7844 

2010 Jun 10 

11734 

LCO/Swope 

WashTl 

050803 


WFC3/IR 

F160W 

906 

2011 Sep 3 

12307 

SwiftfXKY 


050820 

2.612 

ACSAVFCl 

F850LP 

14280 

2006 Jun 8+11 

10551 

HSTIKCSHNVCl 

F850LP 




F625W 

2238 

2006 Jun 5 







F775W 

4404 





050824 

0.83 

WFC3/IR 

F160W 

906 

2011 Jan 18 

12307 

VLT/FORS2 

R 

050904 

6.29 

ACSAVFCl 

F850LP 

4216 

2005 Sep 26 

10616 

HSTmCMOS 

F160W 



NICMOS 

F160W 

15359 

2006 Jul 23 




050908 

3.344 

ACSAVFCl 

F775W 

7892 

2010 Oct 31 

11734 

Gemini-S/GMOS 

r 

051016B 

0.936 

WFC3/IR 

F160W 

906 

2011 Jan 16 

12307 

SwiftrcyoT 

U+B+V 

051022 

0.8 

ACSAVFCl 

F606W 

1560 

2009 Aug 21 

11343 

Chandra 




WFC3/IR 

F160W 

2397 

2009 Oct 12 




060115 

3.53 

ACSAVFCl 

F814W 

7910 

2010 Aug 27 

11734 

VLT/FORSl 

I 

060116 


NICMOS 

F160W 

5120 

2006 Dec 12 

10633 

HSTimCMOS 

F160W 




FI low 

5120 

2006 Dec 11 






ACSAVFCl 

F775W 

4400 

2006 Oct 10 







F850LP 

1650 

2006 Oct 11 




060124 

2.296 

WFC3/IR 

F160W 

1612 

2010 Sep 28 

12307 

SwiftrCNOT 

V 

060206 

4.048 

ACSAVFCl 

F814W 

9886 

2006 Nov 25 

10817 

Palomar 60inch 

I 

060218 

0.033 

WFC3/IR 

F160W 

906 

2010 Oct 12 

12307 

Ff^Sr/ACSAVFCl 

F625W 



ACSAVFCl 

F625W 

2040 

2006 Oct 26 

10551 






F814W 

2040 








F435W 

1530 

2006 Nov 04 







F555W 

1020 





060223 

4.41 

WFC3/IR 

FI low 

8395 

2010 Sep 13 

11734 

SwiftrCNOT 

V 

060418^ 

1.49 

ACSAVFCl 

F625W 

4220 

2006 May 09 

10551 

St\f-HST 





F775W 

4220 








F555W 

4386 

2006 Jul 13 




060502A 

1.51 

WFC3/IR 

F160W 

1209 

2010 Oct 11 

12307 

Gemini-N/GMOS 

r 

060505^ 

0.089 

ACSAVFCl 

F475W 

27774 

2006 May+Jun 

10551 

Gemini-S/GMOS 

g 




F814W 

6840 

2006 Jun 6+7 






WFC3/IR 

F160W 

906 

2011 Aug 3 

12307 



060512^ 

0.443 

WFC3/IR 

F160W 

453 

2011 Feb 23 

12307 

VLT/FORSl 


060522 

5.11 

WFC3/IR 

FI low 

8395 

2010 Oct 17 

11734 

SwiftrCNOT 

White 

060526 

3.21 

ACSAVFCl 

F775W 

7844 

2009 Aug 9 

11734 

VLT/FORSl 

R 

060602A 

0.787 

WFC3/IR 

F160W 

906 

2010 Dec 5 

12307 

SwiftKKV 


060605 

3.78 

ACSAVFCl 

F775W 

7862 

2010 Oct 6 

11734 

SwiftrCNOT 

Wl+V 

060607 

3.082 

ACSAVFCl 

F775W 

7910 

2010 Sep 17 

11734 

Gemini-N/GMOS 

r 

060614 

0.125 

ACSAVFCl 

F606W 

3600 

2006 Sep 8 

10917 

HSTm¥VC2 

F606W 




F435W 

2152 








F814W 

4840 

2006 Oct 31 






WFC3/IR 

F160W 

906 

2010 Oct 8 

12307 



060719 

1.532 

WFC3/IR 

F160W 

1209 

2013 Feb 24 

12949 

VLT/FORS2 

R 




F125W 

1059 





060729 

0.54 

WFC3/IR 

F160W 

906 

2010 Sep 15 

12307 

Gemini-S/GMOS 

i 

060912A 

0.937 

WFC3/IR 

F160W 

906 

2011 Sep 23 

12307 

SwiftrCNOT 

White 

060923A 


WFC3/IR 

F160W 

1198 

2013 Feb 10 

12949 

Gemini-N/NIRI 

K 




F125W 

1348 





060927 

5.6 

WFC3/IR 

FI low 

13992 

2010 Sep 25 

11734 

VLT/FORS2 

I 



NICMOS 

F160W 

10240 

2007 Jun 29 

10926 



061007 

1.261 

WFC3/IR 

F160W 

1209 

2011 Jul 8 

12307 

VLT/FORSl 

R 

061110A 

0.758 

WFC3/IR 

F160W 

906 

2010 Sep 30 

12307 

VLT/FORSl 

R 

061 HOB 

3.44 

ACSAVFCl 

F775W 

7862 

2010 Sep 23 

11734 

VLT/FORSl 

R 

061222A 

2.088 

NICMOS 

F160W 

7680 

2007 Jun 16 

10908 

Gemini-N/NIRI 

K, 
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Table 1 — Continued 


HST Image Afterglow Image 


GRB 

z 

Instrument 

Filter 

Exp. Time 
(s) 

Obs. Date 
(UT) 

Program ID 

Telescope 

Filter 

070125 

1.547 

WFC3/IR 

FI low 

2812 

2010 Apr 24 

11717 

Gemini-N/GMOS 

r 



WFC3/UVIS 

F336W 

2520 





070208 

1.165 

WFC3/IR 

F160W 

1209 

2013 Jan 02 

12949 

Gemini-N/GMOS 

r 




FI low 

1012 





070306 

1.496 

WFC3/IR 

F160W 

1198 

2012 Nov 15 

12949 

VLT/ISAAC 

K. 




F125W 

1348 





070318 

0.836 

WFC3/IR 

F160W 

906 

2010 Dec 31 

12307 

VLT/FORSl 

R 

070508 


WFC3/IR 

F160W 

906 

2011 Feb 22 

12307 

Magellan/MAGIC 

I 

070521 

1.35 

WFC3/IR 

F160W 

906 

2011 Aug 02 

12307 

Swift/XKY 


07072IB 

3.626 

ACSAVFCl 

F775W 

7844 

2010 Nov 13 

11734 

VLT/FORS2 

R 

070802 

2.45 

WFC3/IR 

F160W 

1209 

2013 Feb 22 

12949 

VLT/FORS2 

I 




F105W 

1059 





071010A 

0.98 

WFC3/IR 

F160W 

906 

2010 Oct 29 

12307 



071010B 

0.947 

WFC3/IR 

F160W 

906 

2010 Nov 25 

12307 

Gemini-N/GMOS 

r 

071031 

2.692 

WFC3/IR 

F160W 

1612 

2010 Nov 20 

12307 

VLT/FORS2 

V 

071112C 

0.823 

WFC3/IR 

F160W 

906 

2010 Oct 08 

12307 

Gemini-N/GMOS 

r 

071122 

1.14 

WFC3/IR 

F160W 

1209 

2010 Dec 21 

12307 

Gemini-N/GMOS 

i 

080207 

2.086 

WFC3/IR 

FI low 

2397 

2009 Dec 09 

11343 

Chandra 




WFPC2 

F606W 

1600 

2008 Mar 18 







F814W 

3320 

2009 Mar 20 







F702W 

3600 

2009 Mar 21 






NICMOS 

F160W 

2560 

2008 Apr 05 




080319B‘‘ 

0.937 

WFPC2 

F606W 

3200 

2008 Jul 03 

11513 

Hsir^vvci 





F814W 


2008 Jul 05 




080319C 

1.95 

WFC3/IR 

F160W 

1209 

2010 Sep 19 

12307 

Gemini-N/GMOS 

r 

080325 

1.78 

WFC3/IR 

F160W 

1209 

2012 Dec 31 

12949 

Subaru/MOIRCS 

K. 




F125W 

1359 





080430 

0.767 

WFC3/IR 

F160W 

906 

2011 Jun21 

12307 

SwiftrcyoT 

U+White 

080520 

1.545 

WFC3/IR 

F160W 

1209 

2011 Feb 08 

12307 

SwiftrCNOT 

White 

080603A 

1.688 

WFC3/IR 

F125W 

1397 

2013 Apr 10 

12949 

Gemini-N/GMOS 

r 




F160W 

1348 





080603B 

2.69 

WFC3/IR 

F160W 

1612 

2011 Aug 06 

12307 

SwiftrCNOT 

U+B+V+Wl 

080605'^ 

1.640 

WFC3/IR 

F160W 

2418 

2012+2013 

12307+12949 

VLT/FORS2 

R 




F125W 

1209 

2013 Mar 15 

12949 



080607 

3.036 

WFC3/IR 

F160W 

10791 

2010 Jul 25 

12005 

SwiftrcNOi: 

White+V+W2 

080707 

1.23 

WFC3/IR 

F160W 

1209 

2010 Oct 31 

12307 

VLT/FORSl 

R 

080710 

0.845 

WFC3/IR 

F160W 

906 

2011 Feb 12 

12307 

Gemini-S/GMOS 

r 

080805 

1.505 

WFC3/IR 

F160W 

1209 

2011 Oct 01 

12307 

VLT/FORS2 

I 

080913 

6.7 

WFC3/IR 

F160W 

7818 

2009 Nov 30 

11189 

VLT/FORS2 

z 

080916A 

0.689 

WFC3/IR 

F160W 

906 

2011 Mar 19 

12307 

VLT/FORSl 

R 

080928 

1.692 

WFC3/IR 

F160W 

1209 

2010 Sep 18 

12307 

Gemini-S/GMOS 

r 

081007 

0.530 

WFC3/IR 

F160W 

906 

2010 Nov 30 

12307 

Gemini-S/GMOS 

r 

081008 

1.969 

WFC3/IR 

F160W 

1209 

2011 Sep 04 

12307 

Gemini-S/GMOS 

r 

081109 

0.979 

WFC3/IR 

F160W 

1359 

2013 Mar 09 

12949 






FI low 

1312 





081121 

2.512 

WFC3/IR 

F160W 

1612 

2011 Jan 04 

12307 

SwiftrCNOi: U+B+V+White 

081221 

2.26 

WFC3/IR 

F160W 

1209 

2013 Jun 17 

12949 

Gemini-N/NIRI 

K 




F105W 

1209 





090113 

1.749 

WFC3/IR 

F160W 

2612 

2009 Oct 17 

11840 





ACSAVFCl 

F606W 

2208 

2009 Oct 15 




090404 


WFC3/IR 

F160W 

2612 

2010 Jan 09 

11840 

PdBI 

108 GHz 




F105W 

1359 

2013 Jul 06 

12949 






F125W 

1209 







ACSAVFCl 

F606W 

2208 

2010 Sep 02 

11840 



090407 

1.449 

WFC3/IR 

F160W 

1209 

2010 Sep 15 

11840 

Chandra 




WFC3/UVIS 

F606W 

740 





090417B 

0.345 

WFC3/IR 

F160W 

2612 

2009 Oct 17 

11840 

SwiftlXKT 




ACSAVFCl 

F606W 

1656 

2011 Jan 22 




090418A 

1.608 

WFC3/IR 

F160W 

1209 

2010 Oct 02 

12307 

SwiftrCNOT 

U 

090423^ 

8.2 

WFC3/IR 

F160W 

42997 

2010+2011 

11189 

Gemini-N/NIRI 

H 




F125W 

52117 

2010 




090424 

0.544 

WFC3/IR 

F160W 

906 

2011 May 03 

12307 

Gemini-S/GMOS 

r 

090429Bb 

9.4 

WFC3/IR 

F160W 

4823 

2010 Jan+Feb 

11189 

Gemini-N/NIRI 





F105W 

4823 

2010 Feb 24+28 






ACSAVFCl 

F606W 

2100 

2010 Jan 03 




090618 

0.54 

WFC3/IR 

F160W 

906 

2011 Jan 09 

12307 

Palomar 60inch 

I 

090709A 


WFC3/IR 

F160W 

1359 

2012 Oct 15 

12949 

Subaru/MOIRCS 

K. 




F125W 

1397 





091127 

0.49 

WFC3/IR 

F160W 

906 

2010 Dec 16 

12307 

Gemini-S/GMOS 

i 

091208B 

1.063 

WFC3/IR 

F160W 

1209 

2010 Oct 10 

12307 

Gemini-N/GMOS 

r 

100205A 


WFC3/IR 

F160W 

1209 

2010 Dec 06 

11840 

Gemini-N/NIRI 

K 



WFC3/UVIS 

F606W 

1140 





100316D 

0.059 

WFC3/UVIS 

F336W 

1680 

2011 Aug 03 

12323 

//*Sr/WFC3/UVIS 

F336W 




F225W 

2520 

2011 Aug 02 
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Table 1 — Continued 


GRB 

z 



HST Image 


Afterglow Image 

Instrument 

Filter 

Exp. Time 
(s) 

Obs. Date 
(UT) 

Program ID 

Telescope 

Filter 




F438W 

1230 








F814W 

1260 








F275W 

2520 

2011 Aug 03 







F555W 

1680 








F625W 

720 

2011 Aug 06 







F763M 

1720 







WFC3/IR 

F125W 

1312 

2011 Aug 02 







F160W 

1312 







ACSAVFCl 

FR716N 

2400 

2011 Aug 06 




100413A 

3.9 

WFC3/IR 

F160W 

1209 

2010 Aug 31 

11840 

EVLA 

8.5 GHz 



WFC3/UVIS 

F606W 

752 





100526A 


WFC3/IR 

F125W 

1348 

2013 Jun 30 

12949 

Gemini-N/NIRI 

K 




F160W 

1209 





100615A 

1.398 

WFC3/UVIS 

F606W 

1128 

2010 Dec 16 

11840 

Chandra 




WFC3/IR 

F160W 

1209 





100621A 

0.542 

WFC3/IR 

F160W 

909 

2013 Mar 09 

12949 

VLT/XSHOOTER 

I 




F105W 

748 








F125W 

898 





100905A 


WFC3/IR 

F140W 

5440 

2013 Jul 10 

12247 



110312A 


WFC3/IR 

F160W 

1209 

2011 Nov 18 

12378 

SwiftfXKT 




WFC3/UVIS 

F606W 

IIIO 

2011 Nov 17 




110709B 


WFC3/UVIS 

F606W 

2480 

2011 Nov 08 

12378 

Chandra 




WFC3/IR 

F160W 

2612 

2011 Nov 12 




110731A 

2.83 

WFC3/UVIS 

F625W 

5200 

2012 Aug 22 

12370 

77^rAVFC3/UVIS 

F625W 

120119A 

1.728 

WFC3/IR 

F160W 

1209 

2012 Oct 08 

12949 

Gemini-S/GMOS 

r 




F125W 

1209 






Note. — The table is organized such that the first image listed for each burst is the image for which the corresponding offset and fractional flux measurements 
were used to compile the distribution (Section]^. For bursts with multiple observations we group the observations by instrument. Within an instrument group, 
the observations are listed by date. Only the last epoch of observations are shown for bursts with multiple epochs, unless an earlier epoch was chosen as the 
best host image. Blank entries in the table indicate inheritance of the entry in the previous row. 

^ The afterglow and what we believe to be the host galaxy jPollack et al.|2009t are both detected in the same epoch and offset from each other such that centroid 
measurements are not biased. Therefore the offset (see Section [ tT) can be measured using this single early epoch image. To measure the fractional flux (see 
Section [3^ we used a late epoch image where the afterglow has faded away. 

^ We combined multiple epochs in the same hlter to obtain deep images. For GRB 060505, the afterglow was not detected in early epoch HST observations so 
we combined these observations with a later epoch. 

^ The HST data for GRB 060512 is of poor quality, preventing further analysis. 

^ We combined the F606W and F814W filters to obtain a higher S/N detection of the host. The afterglow image is also a stack of these two hlters at an earlier 
epoch. 
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Table 2 

Offsets, Galaxy Sizes, Magnitudes, and Pec 


GRB 

# Tie Objects 

<^tie 

(") 

<^AG 

(") 

<^host 

(") 

R 

(") 

^phys 

( kpc ) 

Rh 

(") 

Pnorm 

AB Mag ^ 

Pcc 

040701 

4 

0.367 

0.475 








040812 

5 

0.115 

0.046 

0.019 

0.270 ± 0.125 

2.283 ± 1.057 

0.297 

0.908 ± 0.420 

19.450 ± 0.010 

< 0.001 

040916 



0.150 






> 27.1 


040924 

41 

0.005 

0.002 

0.014 

0.286 ± 0.015 

2.197 ± 0.116 

0.188 

1.521 ± 0.080 

23.954 ± 0.024 

0.003 

041006 '’ 

24 

0.004 

< 0.001 

0.007 

0.352 ± 0.008 

2.536 ± 0.061 

0.298 

1.180 ± 0.029 

24.668 ± 0.057 

0.012 

050315 

17 

0.054 

0.003 

0.017 

0.114 ± 0.056 

0.953 ± 0.474 

0.236 

0.481 ± 0.239 

23.808 ± 0.029 

0.006 

050401 

24 

0.044 

0.025 

0.009 

0.082 ± 0.051 

0.638 ± 0.397 

0.170 

0.483 ± 0.300 

25.217 ± 0.077 

0.006 

050406 

7 

0.038 

0.217 

0.006 

0.195 ± 0.220 

1.579 ± 1.788 

0.121 

1.609 ± 1.822 

26.602 ± 0.140 

0.041 

050408 

14 

0.027 

0.002 

0.010 

0.152 ± 0.029 

1.263 ± 0.239 

0.278 

0.545 ± 0.103 

23.700 ± 0.022 

0.008 

050416A 

24 

0.007 

< 0.001 

0.011 

0.037 ± 0.013 

0.260 ± 0.090 

0.160 

0.234 ± 0.081 

23.001 ± 0.009 

0.001 

050525 

74 

0.006 

0.001 

0.002 

0.028 ± 0.007 

0.185 ± 0.046 

0.119 

0.232 ± 0.057 

26.067 ± 0.064 

0.005 

050730 

9 

0.038 

0.014 






> 27.9 


050803 



0.850 








050820 

28 

0.009 

0.003 

0.004 

0.442 ± 0.010 

3.530 ± 0.081 

0.116 

3.808 ± 0.087 

25.196 ± 0.050 

0.010 

050824 

6 

0.036 

0.012 

0.038 

0.473 ± 0.054 

3.597 ± 0.410 

0.402 

1.177 ± 0.134 

23.716 ± 0.042 

0.021 

050904 

9 

0.027 

0.020 






> 27.2 


050908 

15 

0.054 

0.001 

0.003 

0.037 ± 0.055 

0.275 ± 0.405 

0.068 

0.544 ± 0.802 

27.669 ± 0.136 

0.005 

051016B 

4 

0.062 

0.184 

0.013 

0.257 ± 0.194 

2.023 ± 1.532 

0.228 

1.127 ± 0.853 

22.314 ± 0.011 

0.004 

051022 

5 

0.121 

0.015 

0.021 

0.171 ± 0.124 

1.283 ± 0.932 

0.334 

0.512 ± 0.372 

22.020 ± 0.006 

0.002 

060115 

18 

0.045 

0.043 

0.017 

0.280 ± 0.065 

2.042 ± 0.471 

0.125 

2.239 ± 0.517 

27.342 ± 0.137 

0.030 

060116 

10 

0.024 

0.029 






> 27.1 


060124 

7 

0.058 

0.065 

0.025 

0.132 ± 0.091 

1.081 ± 0.745 

0.172 

0.766 ± 0.528 

26.133 ± 0.126 

0.011 

060206 

6 

0.127 

0.021 

0.005 

0.292 ± 0.128 

2.021 ± 0.888 

0.091 

3.211 ± 1.411 

27.647 ± 0.094 

0.035 

060218 

13 

0.012 

< 0.001 

0.044 

0.174 ± 0.045 

0.115 ± 0.030 

0.606 

0.287 ± 0.075 

19.635 ± 0.003 

0.002 

060223 

9 

0.173 

0.041 

0.009 

0.112 ± 0.178 

0.749 ± 1.186 

0.211 

0.532 ± 0.844 

26.534 ± 0.069 

0.026 

060418 



0.003 

0.004 

0.404 ± 0.005 

3.419 ± 0.040 

0.118 

3.425 ± 0.040 

25.855 ± 0.067 

0.019 

060502A 

8 

0.021 

0.009 

0.006 

0.051 ± 0.023 

0.435 ± 0.198 

0.162 

0.317 ± 0.144 

25.864 ± 0.110 

0.007 

060505 

20 

0.029 

0.004 

0.023 

4.307 ± 0.037 

7.160 ± 0.062 

1.195 

3.605 ± 0.031 

19.101 ± 0.001 

0.009 

060522 

4 

0.089 

0.146 






> 28.9 


060526 

23 

0.040 

0.006 






> 27.4 


060602A 



2.550 








060605 

12 

0.115 

0.020 

0.002 

0.183 ± 0.117 

1.304 ± 0.829 

0.109 

1.683 ± 1.069 

27.555 ± 0.148 

0.026 

060607 

33 

0.037 

0.003 






> 28.9 


060614 

21 

0.009 

0.002 

0.011 

0.357 ± 0.014 

0.801 ± 0.031 

0.329 

1.086 ± 0.042 

22.821 ± 0.007 

0.003 

060719 

12 

0.051 

0.062 

0.011 

0.206 ± 0.081 

1.743 ± 0.684 

0.246 

0.837 ± 0.328 

23.400 ± 0.023 

0.006 

060729 

17 

0.039 

0.001 

0.015 

0.338 ± 0.042 

2.149 ± 0.265 

0.308 

1.098 ± 0.135 

23.509 ± 0.028 

0.011 

060912A 

11 

0.148 

0.005 

0.036 

0.654 ± 0.153 

5.155 ± 1.203 

0.713 

0.918 ± 0.214 

21.547 ± 0.013 

0.017 

060923A 

38 

0.039 

0.013 

0.052 

0.291 ± 0.066 

2.460 ± 0.560 

0.696 

0.418 ± 0.095 

23.278 ± 0.038 

0.038 

060927 

18 

0.026 

0.111 






> 28.6 


061007 

5 

0.026 

0.100 

< 0.001 

0.022 ± 0.103 

0.184 ± 0.862 

0.373 

0.059 ± 0.277 

23.621 ± 0.033 

0.013 

061110A 

11 

0.035 

0.014 

0.015 

0.135 ± 0.040 

0.992 ± 0.297 

0.176 

0.765 ± 0.229 

25.114 ± 0.083 

0.007 

061 HOB 

49 

0.025 

0.011 

0.011 

0.061 ± 0.030 

0.449 ± 0.220 

0.104 

0.586 ± 0.288 

27.230 ± 0.105 

0.008 

061222A 

15 

0.026 

0.034 

0.011 

0.044 ± 0.044 

0.368 ± 0.368 

0.193 

0.229 ± 0.229 

24.399 ± 0.050 

0.005 

070125 

5 

0.033 

0.007 






> 28.0 


070208 

6 

0.030 

0.006 

0.025 

0.099 ± 0.040 

0.820 ± 0.327 

0.368 

0.270 ± 0.108 

21.902 ± 0.009 

0.005 

070306 

14 

0.043 

0.005 

0.011 

0.090 ± 0.045 

0.760 ± 0.380 

0.209 

0.430 ± 0.215 

21.860 ± 0.007 

0.002 

070318 

8 

0.057 

0.007 

0.019 

0.109 ± 0.060 

0.827 ± 0.458 

0.205 

0.529 ± 0.293 

24.100 ± 0.034 

0.005 

070508 

16 

0.025 

0.004 

0.020 

0.423 ± 0.033 

3.577 ± 0.277 

0.592 

0.714 ± 0.055 

22.871 ± 0.031 

0.027 

070521 



0.850 








070721B 

17 

0.042 

0.032 

0.006 

0.045 ± 0.054 

0.327 ± 0.387 

0.057 

0.793 ± 0.939 

28.442 ± 0.165 

0.011 

070802 

10 

0.039 

0.007 

0.027 

0.132 ± 0.048 

1.069 ± 0.390 

0.348 

0.379 ± 0.138 

23.994 ± 0.043 

0.014 

071010A 

11 

0.127 

0.300 








071010B 

7 

0.022 

0.001 

0.010 

0.109 ± 0.025 

0.860 ± 0.195 

0.257 

0.424 ± 0.096 

22.635 ± 0.016 

0.004 

071031 

7 

0.068 

< 0.001 






> 27.4 


071112C 

10 

0.020 

0.007 

0.017 

0.202 ± 0.027 

1.533 ± 0.206 

0.398 

0.508 ± 0.068 

23.674 ± 0.046 

0.016 

071122 

10 

0.032 

0.033 

0.025 

0.075 ± 0.053 

0.616 ± 0.433 

0.393 

0.191 ± 0.134 

22.678 ± 0.017 

0.009 

080207 

2 

0.381 

0.206 

0.066 

0.769 ± 0.438 

6.402 ± 3.650 

0.683 

1.126 ± 0.642 

23.225 ± 0.025 

0.047 

080319B 

18 

0.011 

0.001 

0.117 

0.215 ± 0.118 

1.697 ± 0.929 

0.147 

1.465 ± 0.802 

26.679 ± 0.200 

0.015 

080319C 

10 

0.018 

0.008 

0.046 

0.855 ± 0.050 

7.177 ± 0.417 

0.638 

1.341 ± 0.078 

22.018 ± 0.011 

0.023 

080325 

60 

0.022 

0.041 

0.059 

0.697 ± 0.075 

5.890 ± 0.635 

0.443 

1.574 ± 0.170 

22.537 ± 0.016 

0.017 

080430 

4 

0.035 

0.015 

0.030 

0.108 ± 0.049 

0.800 ± 0.361 

0.282 

0.384 ± 0.173 

24.595 ± 0.071 

0.011 

080520 

17 

0.106 

0.146 

0.024 

0.468 ± 0.182 

3.966 ± 1.543 

0.223 

2.100 ± 0.817 

23.965 ± 0.032 

0.011 

080603A 

22 

0.027 

< 0.001 

0.013 

0.089 ± 0.030 

0.753 ± 0.255 

0.202 

0.440 ± 0.149 

22.884 ± 0.011 

0.003 

080603B 

6 

0.117 

0.034 






> 26.7 


080605 

38 

0.029 

0.009 

0.008 

0.081 ± 0.031 

0.682 ± 0.266 

0.198 

0.407 ± 0.159 

22.467 ± 0.008 

0.002 

080607 

5 

0.230 

0.056 

0.024 

0.681 ± 0.238 

5.229 ± 1.827 

0.419 

1.626 ± 0.568 

24.372 ± 0.027 

0.036 

080707 

6 

0.047 

0.024 

0.009 

0.089 ± 0.054 

0.739 ± 0.449 

0.263 

0.338 ± 0.205 

22.950 ± 0.017 

0.005 

080710 

6 

0.035 

< 0.001 






> 26.1 


080805 

15 

0.038 

0.006 

0.026 

0.466 ± 0.047 

3.944 ± 0.396 

0.275 

1.695 ± 0.170 

23.249 ± 0.020 

0.010 

080913 

13 

0.035 

0.044 






> 28.0 


080916A 

12 

0.056 

0.022 

0.010 

0.011 ± 0.061 

0.075 ± 0.432 

0.240 

0.044 ± 0.254 

22.762 ± 0.016 

0.004 

080928 

16 

0.021 

0.012 

0.016 

1.705 ± 0.029 

14.437 ± 0.245 

0.442 

3.858 ± 0.066 

21.990 ± 0.010 

0.035 

081007 

11 

0.023 

< 0.001 

0.017 

0.111 ± 0.029 

0.700 ± 0.182 

0.277 

0.402 ± 0.105 

24.489 ± 0.059 

0.011 

081008 

28 

0.021 

0.002 

0.018 

1.658 ± 0.028 

13.902 ± 0.236 

0.236 

7.027 ± 0.119 

23.591 ± 0.028 

0.066 
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Table 2 — Continued 


GRB 

# Tie Objects 

^tie 

(") 

(") 

<^host 

(") 

R 

(") 

^phys 

(kpc) 

Rh 

(") 

Rnorm 

AB Mag^ 

Tec 

081109 



0.850 










081121 

6 

0.12 

0.014 

0.011 

0.262 


0.122 

2.112 ±0.981 

0.198 

1.323 ±0.615 

24.923 ± 0.054 

0.010 

081221 

11 

0.047 

0.016 

0.055 

0.388 


0.074 

3.195 ±0.609 

0.409 

0.949 ±0.181 

23.156 ±0.025 

0.013 

090113 



0.910 










090404 

4 

0.115 

0.361 

0.044 

0.556 


0.381 

4.709 ± 3.229 

0.690 

0.806 ± 0.553 

23.370 ±0.018 

0.045 

090407 

5 

0.134 

0.387 

0.021 

0.172 


0.410 

1.451 ±3.463 

0.336 

0.511 ± 1.220 

22.953 ± 0.019 

0.026 

090417B 



0.850 










090418A 

5 

0.209 

0.014 

0.017 

0.096 


0.211 

0.812 ± 1.785 

0.214 

0.448 ± 0.984 

23.687 ± 0.025 

0.009 

090423 

6 

0.031 

0.005 








>28.2 


090424 

9 

0.021 

0.005 

0.030 

0.407 


0.037 

2.594 ± 0.235 

0.407 

1.000 ±0.090 

21.260 ±0.007 

0.005 

090429B 

10 

0.053 

0.012 








>28.1 


090618 

12 

0.038 

0.001 

0.025 

0.693 


0.046 

4.402 ± 0.292 

0.519 

1.335 ±0.089 

22.464 ± 0.020 

0.020 

090709A 

25 

0.028 

0.021 

0.007 

0.283 


0.036 

2.393 ± 0.305 

0.256 

1.105 ±0.141 

24.213 ± 0.038 

0.010 

091127 

16 

0.040 

0.005 

0.028 

0.221 


0.049 

1.335 ±0.296 

0.411 

0.538 ±0.119 

22.700 ± 0.021 

0.011 

091208B 

7 

0.019 

0.001 








>26.7 


100205A 

12 

0.034 

0.028 








>26.7 


100316D" 

15 

0.007 

<0.001 

<0.001 









100413A 

8 

0.098 

0.600 










100526A 

10 

0.047 

0.008 

0.008 

0.194 


0.049 

1.645 ±0.411 

0.188 

1.034 ±0.258 

24.766 ± 0.033 

0.007 

100615A 

3 

0.318 

0.028 

0.004 

0.246 


0.320 

2.073 ± 2.694 

0.083 

2.963 ± 3.850 

25.092 ± 0.033 

0.034 

100621A 

6 

0.062 

0.025 

0.015 

0.045 


0.069 

0.284 ± 0.439 

0.262 

0.17 ± 0.263 

21.316 ±0.006 

0.002 

100905A 



0.910 










110312A 



0.850 










110709B 

1 

0.318 

0.057 

0.005 

0.067 


0.323 

0.564 ± 2.737 

0.073 

0.912 ± 4.430 

26.714 ± 0.080 

0.111 

110731A 

211 

0.006 

0.001 

0.008 

0.182 


0.010 

1.423 ±0.082 

0.121 

1.502 ±0.087 

25.236 ± 0.043 

0.005 

120119A 

49 

0.013 

<0.001 

0.012 

0.012 


0.017 

0.099 ±0.146 

0.202 

0.058 ± 0.086 

23.480 ± 0.022 

0.003 


Note. — For bursts with multiple HST images of the host galaxy, the image used to make the measurements listed here corresponds to the first image listed 
for such bursts in Table|T] We show the 3cr upper limits for bursts with no host detection. Bursts with no listed offset or upper limit correspond to those bursts 
for which we were not ^e to obtain a suitable afterglow image. 

^ Magnitudes are corrected for Galactic extinction jSchlafly & Finkbeiner|2Ql l} 

^ The SN associated with GRB 041006 is detected in the final epoch. The SN and host center are sufficiently offset from each other such that centroid 
measurements are not biased. However, measurement of fractional flux is not possible. 

^ GRB 100316D occurred in a low-redshift galaxy with irregular morphology. The complexity of the host system prevents a meaningful definition of the host 
center. We therefore do not make an offset measurement for this burst. 
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Table 3 

Fractional Flux and Ratio of Error Circle to Galaxy Area 


GRB 

Fractional Flux 

Ratio 

GRB 

Fractional Flux 

Ratio 

040701 



070802 

0.692 

1.30 X 10-2 

040812 

0.423 

1.74 X 10"' 

071010A 



040916 



071010B 

0.825 

7.34 X 10-2 

040924 

0.783 

8.21 X 10-^ 

071031 



041006 


1.80 X 10-^ 

071112C 

0.774 

2.83 X 10-2 

050315 

0.693 

5.25 X 10-2 

071122 

0.951 

1.37 X 10-2 

050401 

0.799 

8.86 X 10-2 

080207 

0.207 

4.02 X 10-' 

050406 

0.110 

3.31 

080319B 

0.598 

5.65 X 10-2 

050408 

0.724 

9.48 X 10-2 

080319C 

0.599 

9.53 X 10-+ 

050416A 

0.896 

1.91 X 10-2 

080325 

0.134 

1.10 X 10-2 

050525 

0.914 

2.61 X 10-2 

080430 

0.645 

1.82 X 10-2 

050730 



080520 

0.117 

6.55 X 10-' 

050803 



080603A 

0.880 

1.79 X 10-2 

050820 

0.606 

6.69 X 10-2 

080603B 



050824 

0.324 

8.91 X 10-2 

080605 

0.820 

2.35 X 10-2 

050904 



080607 

0.177 

3.19 X 10-' 

050908 

0.316 

6.31 X 10-‘ 

080707 

0.836 

4.03 X 10-2 

051016B 

0.329 

7.25 X 10-' 

080710 



051022 

0.556 

1.33 X 10-' 

080805 

0.517 

1.96 X 10-2 

060115 

0.081 

2.48 X 10-' 

080913 



060116 



080916A 

0.930 

6.28 X 10-2 

060124 

0.427 

2.57 X 10-' 

080928 

0.009 

2.99 X 10-2 

060206 

0.032 

2.00 

081007 

0.792 

6.89 X 10-2 

060218 

0.872 

3.92 X 10-^ 

081008 

0.000 

7.99 X 10-2 

060223 

0.292 

7.10 X 10-' 

081109 



060418 

0.000 

6.46 X 10-^ 

081121 

0.197 

3.72 X 10-' 

060502A 

0.689 

1.99 X 10-2 

081221 

0.761 

1.47 X 10-2 

060505 

0.987 

6.00 X 10-^ 

090113 



060522 



090404 

0.183 

3.02 X 10-' 

060526 



090407 

0.245 

1.49 

060602A 



090417B 



060605 

0.103 

1.15 

090418A 

0.367 

9.58 X 10-' 

060607 



090423 



060614 

0.459 

7.85 X 10-^ 

090424 

0.740 

2.81 X 10-2 

060719 

0.581 

1.07 X 10-' 

090429B 



060729 

0.246 

1.60 X 10-2 

090618 

0.260 

5.36 X 10-2 

060912A 

0.525 

4.31 X 10-2 

090709A 

0.461 

1.87 X 10-2 

060923A 

0.980 

3.49 X 10-2 

091127 

0.833 

9.62 X 10-2 

060927 



091208B 



061007 

0.778 

7.67 X 10-2 

100205A 



061110A 

0.419 

4.59 X 10-2 

100316D 

0.879 


061 HOB 

0.393 

6.90 X 10-2 

100413A 



061222A 

0.955 

4.92 X 10-2 

100526A 

0.531 

6.43 X 10-2 

070125 



100615A 

0.051 

14.8 

070208 

0.928 

6.91 X 10-2 

100621A 

0.914 

6.51 X 10-2 

070306 

0.798 

4.29 X 10-2 

100905A 



070318 

0.781 

7.85 X 10-2 

110312A 



070508 

0.517 

1.83 X 10-2 

110709B 

0.055 

19.6 

070521 



110731A 

0.029 

2.53 X 10-2 

070721B 

0.211 

8.58 X 10-' 

120119A 

1.000 

4.14 X 10-2 
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Table 4 

Summary of KS Test p-values 


LGRBs - This Work 

Exp. Disk 

Type Ib/c SNe^ 

Type Ic SNe 

Type II SNe 

CCSNe 

SLSNe 

LGRBs - F06 LGRBs - B02 

Physical Offsets 


1.33x10“'’ 


1.53x10“’ 



0.84 

Host-Normalized Offsets 

2.6x10“^ 

1.3x10“2 


5.2x10“5 


0.17 

0.21 

Fractional Fluxes 



0.79 

5.4x10“5 

3.0x10“’ 

0.52 

0.09 


Note. — The values reported here for the physical offsets, host-normalized offsets, and fractional fluxes include only bursts with ^1-0 kpc, ^ 

0.5, and error circle to galaxy area ratio <0.1, respectively. Comparison samples of Type Ib/c and II SNe (physical offsets). Type Ib/c and II SNe (host- 
norma l ized offsets) . Type Ic and II SNe (fraction a l fluxes), CCSNe, SLSNe, LGRBs - F06, a nd L GRBs - B02 are from |Prieto et al.M2008),|Kelly & Kirshnerj 
j2012) , [Kelly et al.|j2008) , [^ensson et ah] ( |2010) , |Lunnan et al.H2015) , |Fruchter et al.lpOObf and |Bloom et al.jpQO^ , respectively. 

^ For the host-normalized offsets, 80% of LGRB Monte Carlo synthetic distributions have p-values > 0.05 when compared to Type Ib/c SNe. 


Table 5 

Summary of LGRB Offset and Fractional Flux Medians 



Sample Median 

Mean of Monte Carlo Medians 

90% Confldence Interval 

Physical Offsets 

1.27 kpc 

1.34 kpc 

1.17-1.51 kpc 

Host-Normalized Offsets 

0.54 

0.67 

0.58 - 0.76 

Fractional Fluxes 

0.75 




Note. — The values reported here for the physical offsets, host-normalized offsets, and fractional fluxes include only bursts with ^1-0 kpc, < 

0.5, and error circle to galaxy area ratio <0.1, respectively. 




























